Effective action for vortex dynamics in clean d-wave superconductors 
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We describe the influence of the gapless, nodal, fermionic quasiparticles of a two-dimensional d- 
wave superconductor on the motion of vortices. A continuum, functional formalism is used to obtain 
the effective vortex action, after the fermions have been integrated out. At zero temperture (T), 
the leading terms in the vortex action retain their original form, with only a finite renormalization 
of the vortex effective mass from the fermions. A universal "sub-Ohmic" damping of the vortex 
motion is also found. At T > 0, we find a Bardeen-Stephen viscous drag term, with a universal 
co-efficient which vanishes as ~ T 2 . We present a simple scaling interpretation of our results, in 
which quantum-critical Dirac fermions respond to a moving point singularity. Our results appear 
to differ from those of the semiclassical theory, which obtains more singular corrections to a vortex 
mass appearing in transport equations. 
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I. INTRODUCTION 

The cuprate superconductors behave like conventional 
BCS superconductors in many respects, albeit with a d- 
wave symmetry of the Cooper pair wavefunction. The 
BCS wavefunction provides a reasonable description of 
the ground state. The elementary excitations above the 
ground state are 

(i) fermionic, S — 1/2, Bogoliubov quasiparticles, 

(ii) plasmons (or "phase fluctuations"), and 
(Hi) vortices with magnetic flux hc/(2e). 

For the traditional "low temperature" superconduc- 
tors, the vortices are usually treated as a classical, con- 
figuration of the superconducting order parameter, which 
forms a background for the quantum fluctuations of the 
fermionic quasiparticles and the Cooper pairs: the quan- 
tum zero-point motion of the vortices themselves is neg- 
ligibly small. 

Quantum fluctuations of the vortices are expected to 
play a more fundamental role in the cuprate super- 
conductors, requiring that the vortices be treated as 
bona fide quasiparticles and elementary excitations in 
their own right. There are a variety of experimental 
indicationsA*2iiMi£i£il that the superconducting state is 
proximate to an insulating ground state. Vortices and 
anti-vortices are expected to proliferate in the vicinity of 
such a transition, and so the energy gap towards creation 
of vortex/anti- vortex pairs in the superconductor be- 
comes vanishingly small as the superconductor-insulator 
transition is approached. Moreover, the superconductor 
coherence length, £, is at most a few lattice spacings, 
which suggests that the effective mass of the vortices is 
small. 

It is the purpose of this paper to discuss some features 
of the effective action controlling interplay between the 
three classes of elementary excitations (noted above) of 
a two-dimensional d-wave superconductor. If we neglect 
the fermionic S = 1/2 quasiparticles, the required action 
can be readily deduced from the theory of boson-vortex 
duality&Siifl, as we will review below. We will primarily 
be concerned here on the influence of the gapless, nodal, 
fermionic quasiparticles on the vortices. The interaction 




FIG. 1: Scales of a vortex. The core size is of the order 
of superconducting coherence length £, while the extent of 
supercurrent circulations is roughly equal to the much larger 
London penetration depth A. This paper will work in the 
limit in which A — > oo. 



between the fermions and the vortices is non-local, with 
no natural expansion parameter, and so its description 
poses a strong-coupling problem. Here, we will consider 
a single, isolated, slowly moving vortex, and exactly eval- 
uate the functional integral over the fermionic degrees 
of freedom by a combination of analytic and numerical 
techniques; this yields a renormalization of the effective 
action for the vortices and the plasmons. Our primary 
result will be that despite the absence of an energy gap in 
the fermionic spectrum, the influence of the fermionic de- 
grees of freedom is relatively innocuous, and at zero tem- 
perature (T) the leading terms in the vortex action retain 
the form deduced while ignoring the fermionic quasipar- 
ticles. 

The interplay between fermions and moving vortices 
has also been examined previously in semiclassical frame- 
work by a number of workers These theories are 
expansion in the ratio of the superconducting pairing en- 
ergy to the Fermi energy, a ratio which is not small in 
the cuprate superconductors. We will not make any as- 
sumptions on the value of this ratio in the present paper. 
The results of semiclassical theory are expressed in a dif- 
ferent framework of transport equations, and this makes 
direct comparison with our results difficult. An impor- 
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tant characteristic of the semiclassical results is that the 
effective mass of a vortex (as defined in their transport 
equations) in a d-wave superconductor diverges linearly 
with the spacing between vortices. We do not find any 
corresponding divergence in the vortex mass in our func- 
tional formalism. The reasons for this potential discrep- 
ancy between our and the semiclassical results are not 
completely clear to us. We speculate it may be due to 
the breakdown of the semiclassical approximation in de- 
scribing the low energy nodal fermionic quasiparticles. 
Such fermionic excitations are not semiclassical, but, 
with their scale-invariant Dirac spectrum, better char- 
acterized as a quantum-critical system^ 

We begin our presentation by reviewing the action for 
vortices and plasmons as derived in theories of the boson- 
vortex duality?Si2iifl while ignoring fermionic excitations 
of the superconductor. We use a first-quantized approach 
for the vortex degrees of freedom, and treat each vortex 
as a quantum particle at spacetime position r v (r), where 
r is imaginary time; the quantum mechanics of the vor- 
tices is described by a functional integral over r v (r). The 
density or "phase" fluctuations of the superconductor are 
represented by a U(l) gauge field A^(r,r), where the 
spacetime index p extends over x, y, r: the 'magnetic' 
flux of this gauge field is proportional to the particle den- 
sity of the superconductor, while the 'electric' field deter- 
mines the local superflow velocity. A vortex (anti- vortex) 
is minimally coupled to the U(l) gauge fields with charge 
+1 (— 1). So a single vortex and the phase fluctuations 
are described by the action S = S v + Sa where 



S v — 



(It 



m° ( dr v (r) 
2 \ dr 



,ir T (r) jh, , ss 



+ 1 



J d 2 rdTA^{r,T)J v ^{r,T) 



(1) 



is the vortex action and its coupling to the U(l) gauge 
field. Here m° is the bare vortex mass (which can 
be viewed as arising from high energy degrees of free- 
dom which have been integrated out), A represents the 
'background' flux which is responsible for the Magnus 
force on the vortices. In a Galilean-invariant superfluid, 
V x A° = z27r(density of Cooper pairs). In the pres- 
ence of a lattice, when the superconductor is proximate 
to a Mott insulator, it has been arguediSiiLi2iiS that the 
Magnus force is reduced, and the density of pairs in the 
Mott insulator has to be subtracted from the density of 
Cooper pairs. In the last terms in S, J v ^ is the vortex 
3-current, and its temporal and spatial components are 
given by 



J VT (r,r) = 5(r-r v (r)) 
J„(r,r) = S{r-r v (r)) 



dr v (r) 
dr 



(2) 



The phase fluctuations of the superconductor are con- 
trolled by the action Sa which, in the Coulomb gauge 



V • A = 0, is 

, drkdu I 1 



8tt 3 \8tt 2 p s 



[k 2 \A T (k,u)\ 2 + LU 2 \A l {k,U J ){< 



4tt 



kikj 
5 - - I? 



Ai{-k,-w)Aj{k,uj) . (3) 



We have Fourier transformed to momentum, k, and imag- 
inary frequencies, u>, e* = 2e is the charge of a Cooper 
pair, and p s is the superfluid stiffness in units of en- 
ergy (the London penetration depth, A, is related to 
p s by p s — h 2 c 2 d/(Wire 2 A 2 ), where d is the interlayer 
spacing). The first term in Sa is responsible for the 
log interaction between the vortices; the coefficient of 
this term is fixed by requiring that the log interaction 
have the required prefactor. The second term represents 
the electrical Coulomb interaction between the charged 
Cooper pairs; its coefficient is fixed by the knowledge that 
(V x A)/(2n) is the fluctuation in the number density of 
pairs of charge e*. The 'photon' excitations of the full 
action Sa are the plasmons of the two-dimensional super- 
conductor with dispersion lo ~ y/k, and these plasmons 
can be 'radiated' by a moving vortex. 

It is worth pointing out here a noteworthy aspect of 
the action for vortices and plasmons presented here in 
Eqs. JU and 10. Although obtained by an approximate 
and uncontrolled boson- vortex duality mapping, most of 
the coupling constants in the final action can be fixed 
rather precisely. The primary unknown is the vortex 
mass, m v , and we will describe microscopic contributions 
to its value this paper. 

Clearly, the physics described by S v + Sa is similar to 
the quantum electrodynamics of a charged particle. The 
coupling to the phase fluctuations will renormalize the 
vortex mass, m v , and lead to 'radiation damping' in the 
motion of a vortex. We review this physics in Section ITU 
For a neutral superfluid, it has been arg ucd 20 - 21 that 
the phase fluctuations lead to a logarithmically divergent 
"dynamic mass" of the vortex; we show how the analog 
of this result obtains in our formalism in Section [TJ] Sec- 
tion [Tj] also considers the case of a charged superfluid, 
with plasmon excitations as in Eq. for this case, we 
find that the phase fluctuations are less singular and lead 
only to a finite renormalization of the vortex "dynamic 
mass" . 

Our purpose here is to extend S v + Sa to include 
the remaining elementary excitations of the d-wave su- 
perconductor, the nodal fermionic quasiparticles. (The 
complete action controlling the interactions between all 
the elementary excitations of a d-wave superconductor 
^22,23,24,25 § V +S A + S^ + S cs specified in Eqs. (HIJ>, ©, 
(|63|l . and (|u4f> ). We will then integrate out these quasi- 
particles, and associated gauge fields (including A^), in 
the presence of a single, slowly moving vortex, and ob- 
tain resulting renormalization of S v . For small r v (w), the 
effective vortex action can be written in the form (apart 
from terms responsible for the Magnus force associated 
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with A° in Eq. ©) 

s v = J ^F||MKH| 2 + ... (4) 

For small to, we find for a charged d-wave superconductor 
at T = that the contribution of the nodal quasiparticles 
to F\\ (u>) can be written as 

v 2 F F\\(u) = u 2 AFx(y±/v F ) - \u>\ 3 F 2 (v A /v F ), (5) 

where v F and i>a are the two velocities of the nodal quasi- 
particles, A is a high energy cutoff, and F\ t 2 are functions 
of the ratio of the velocities. Comparing with Eq. Q), 
we see that the first term can be interpreted as finite 
renormalized mass m v = AFi(va/v f )/(2v f ) of the vor- 
tex. The function Fi is not universal, and should clearly 
depend upon the nature of the cutoff; however, no other 
coupling constant appears in the value of F\ . The second 
| a; | 3 term in Eq. © represents a "sub-Ohmic" damping 
of the vortex motion. We will show that this damping 
has a universal strength i.e. F 2 is a universal function of 
va/v f , and so the damping is entirely determined by the 
two quasiparticle velocities. (For completeness, we note 
that in a neutral <i-wave superfluid, as we will demon- 
strate in Section m there is an additional w 2 ln(A/|w|) 
term in -Fji(w): this term is the manifestation, in the 
present formalism, of the logarithmically divergent "dy- 
namic mass" of the vorte»22*2i, noted above.) 

It is interesting to note here that the expression for 
Fn (lu) in Eq. (J5J could have been deduced very simply 
and directly by dimensional arguments, after assuming 
that the nodal fermions are properly thought of as a 
'quantum critical' system with dynamic exponent z — 1, 
which then responds to the moving vortex. In this scaling 
argument, the vortex position r v (r) scales as a length, 
and so dim[r v (r)] = — 1. After a Fourier transformation, 
dim[r v (o;)] = —2. Counting dimensions in the action, we 
conclude then that dimfi^ii (a>)] = 3. Also, we know from 
translational invariance that F\\ (0) = 0. So apart from a 
universal |w| 3 term, F\\(lu) can have a cutoff-dependent, 
analytic to 2 term. The dependence of Fy (lu) on the veloc- 
ities follows from ensuring consistency with engineering 
dimensions. The success of this argument demonstrates 
the importance of thinking about the fermions in terms 
of the scale-invariant critical field theory that obeys hy- 
perscaling properties. As we noted earlier, we believe 
that it is this criticality which is not properly captured 
by the semiclassical analysis 

This scaling argument, and our explicit computations, 
show that no \lu\ "Ohmic" damping term is generated 
in (lu) at T = 0. If present, such a term would con- 
tribute the familiar Bardeen-Stephen dampingSi to the 
vortex equations of motion. If the vortex mass had the 
degree of infrared divergence claimed in the semiclassical 
theoryii*i2iiiii^ i.e. a mass which diverges linearly with 
length scale, then this would be equivalent, in our lan- 
guage, to a scaling dimension dim^j (lu)] = 1, and would 
imply that a \lo\ term should have appeared in our vortex 



action. However, because we find dim[Fy (lu)] — 3, all our 
integrals have a degree of convergence in the infrared to 
allow safe expansion in powers of lu upto order lu 2 , and 
to prevent the appearance of a \lu\ term. 

In the real experimental system, our results show that 
any Bardeen-Stephen damping must come from other 
sources: a non-zero T, or from impurities which induce 
a non-zero density of fermionic states at zero energy. In 
particular, we show in Section IV Dl that at T > 0, the 
function Fn (lu) contains a Bardeen-Stephen term 

F||M = |M + ... (6) 

where the viscous drag co-efficient is universal in the 
clean limit: 

V=y^Fv(v A /v F ), (7) 

V F 

with F v a universal function of the velocity ratio. We 
are able to compute one contribution to rj exactly: 
that due to the statistical phase factor acquired by the 
fermions upon encircling the vortex — the result appears 
in Eq. 1)58(1 . There is a second contribution due to the 
Dopplcr shift which requires a numerical analysis. 

We will begin by setting up the general 
formalism 27 - 28 ' 29,30 under which the effective action 
of a single vortex can be derived in Section ITTll The na- 
ture of the coupling between the fermionic quasiparticles 
and the vortices and the plasmons will then be reviewed 
in Section IIVI here we will follow the framework of 
Tesanovic and othersi 31 i 32 i 33 i 34 i 35 

With the basic formalism in place, it then remains 
to explicitly evaluate the functional determinant of the 
nodal fermions. There is no natural small parameter in 
this determinant, and so an exact evaluation is required. 
Nevertheless, it is useful to examine the structure of the 
perturbation theory order-by-order in the non-linearities; 
this will be done in Section lVll We will find that the per- 
turbation theory has the scaling structure discussed near 
Eq. Q. 

The non-perturbative evaluation of the fermion deter- 
minant is presented in Sections IVl and IVTT1 In Section [V] 
we neglect the 'Doppler shift' of the fermion energy in 
the presence of a superflow, and retain only the 'Berry 
phase' acquired by the fermions upon encircling a vortex; 
we show that an analytic evaluation of the effective mass 
renormalization is possible in this limit. Section IVIII ac- 
counts for the Doppler shift in a numerical computation. 
These non-perturbative results are found to be entirely 
consistent with the expectation of the perturbation the- 
ory discussed in Section IVII in that they also has the 
scaling structure of Eq. (JSJ . 

Finally, conclusions and implications of our results ap- 
pear in Section fVIIII and some technical details and sup- 
porting calculations appear in the appendices. 
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II. "DYNAMIC" VORTEX MASS 

This section takes a slight detour. We will review ear- 
lier computations on the influence of phase fluctuations 
on the vortex mass. Consideration of the influence of 
nodal quasiparticles will begin in Section Ifffl 

Popov 20 and DuanSi have argued that vortices acquire 
a "dynamic" mass from their coupling to quantum fluc- 
tuations of the phase of the superconducting order, and 
that this contribution diverges logarithmically for a neu- 
tral superfluid. We now present the analog of this effect 
in our formalism, and show further that the divergence 
is absent in charged superfLuids. 

We begin with a discussion of neutral superfluids. For 
this case, in the action for the phase fluctuations in 
Eq. J2J), the pre- factor k in the last term is replaced by k 2 : 
this ensures only short-range interaction proportional to 
the square of the local density fluctuation (V x A)/(2ir). 
After a suitable rescaling of coordinates to absorb factors 
of the velocity of spin- waves, the resulting action can be 
written in the compact 'relativistic' form 

c f cftp 

$An = 2 8tt3^ (<W 2 - VvVv) A v , (8) 

where c is a coupling constant and — (uj, k) is a space- 
time 3-momentum. Now we can use S v in Eq. and 
integrate out the A^ fluctuations under the action S An . 
In this manner, we obtain the following additional term 
in the vortex effective action 



SvA = IJdrdr'J 



d 2 kdw e -'M T - T ') 
~ (w 2 + k 2 ) 



1 



dr v (r) dr v (T r ) 



dr 



dr' 



(9) 



Integrating over k and u>, we obtain for large r — t' 

1 



T7? ( 10 ) 

[(r v (r)-r v (r')) 2 + (T-T') 2 ] 1/2 

As we will repeatedly do in all the computations of this 
paper, we expand this effective action for r(r) in powers 
of a slowly varying r v (w). This yields one term with the 
structure 



s vA ~ J cWKHl 2 jH dr^ 



(ii) 



(where A is a high energy cutoff), and another with the 
structure 



SvA ~ / duj\r v (uj 



„ , (12) 

A -l T J 



Both terms yield a contribution to Fy (uj) in Eq. (0} which 
is ~ lu 2 ln(A/|w|). This logarithm is an indication of 



the d iverge nce of the dynamic mass, as found earlier in 
Refs. l2(l2lL 

Let us now extend the above analysis to the charged 
case of interest to the cuprates, where the A^ action is 
in Eq. @. Now it is convenient to work in the Coulomb 
gauge, when the only potentially singular contribution 
to the vortex action comes from the fluctuations of the 
spatial components A. In this case, Eq. © is replaced 
by 



S, 



vA 



e" 

dr v i(r) 
dr 



drdr' 



d 2 kduj 



-iuj{r—r') 



8tt 3 (k + uj 2 /(2irp s e* 2 )) 

hkj\ dr VJ (T r ) ik-(r v (r)-r v (r')) 
k 2 / dr' 



(13) 

As above, we expand in powers of r v (u>), and perform the 
momentum and frequency space integrals to obtain an ex- 
pression like that in Eq. (H), but with the 1/r under the 
r integral replaced by 1/r 3 . This r integral is convergent 
in the infrared, and hence we obtain a uj 2 term in F\\ (uj) 
with a finite co-efficient i.e. a finite renormalization of 
the vortex mass. 



III. EFFECTIVE VORTEX ACTION 

We begin by deriving a general effective action for a 
single vortexi2L2S422i2£ We will focus on the contributions 
of the fermionic quasiparticles, with the aim of renor- 
malizing the action S v + S A presented in Section [I] The 
vortex will be treated as a singular configuration of ex- 
ternal gauge and supercurrent fields in the quasiparticle 
Hamiltonian. We write a path-integral for quasiparticles 
and immediately integrate them out, which gives us the 
vortex action as a functional of the time-dependent vor- 
tex location. This is done in an expansion in the vortex 
velocity. Formally, for this procedure to be well defined, 
quasiparticles must be massive. Therefore, we will ini- 
tially give a finite mass to the nodal quasiparticles, and 
send it to zero later when it becomes safe. 

The vortex action is the logarithm of the fermion de- 
terminant in the path-integral: 



S v 



~tr logG 



(14) 



Here, the time-ordered quasiparticle Green's function in 
imaginary time satisfies: 



^ + H) G(r,T;r',T') = 5{t - r')6 2 {r - r')J . ! 1.5) 



where / is the unit-matrix for spin degrees of freedom. In 
a d-wave superconductor there are four flavors of nodal 
quasiparticles - their contributions must be calculated 
separately and added to the vortex action. 

The quasiparticle Hamiltonian H, which we discuss in 
the next section, depends on the vortex position r v (r). 
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We expand it to the second order in small vortex dis- 
placements from the origin: 



H = H Q + (r v V v )H Q + ^(r v V v ) 2 H { 
= H Q + AH(r v ) , 



(16) 



where V v = -^r- We also expand the Green's function, 
and write in the operator notation: 

AH, 



G- 1 = Gq 1 



(17) 



with Go being the unperturbed Green's function (vortex 
sitting at the origin). Now expand (|14l) : 

5 V = -trlogG^ 1 - tr(GoAH) + ^tr(G AH) 2 + ■■■ . 

(18) 

The first term is a constant and we drop it. In order to 
evaluate the other terms, we substitute the solution of 
(fH))l for the unperturbed Green's function: 



G (r,r;r',r' 



— ilo{t— t ) 



>„(r)^(r') 



(19) 

The wavefunctions ip n and energies e n must be solutions 
of the Schrodinger equation in position representation: 



H (r)ip n (r) = e n ip n (r) 



(20) 



and frequencies u> take discrete values (2k + l)TrT at finite 
temperature T = 1//3. Clearly, the vortex dynamics can 
be derived entirely from the knowledge of quasiparticle 
energies and wavefunctions in presence of the vortex at 
rest. 

The second term of (IT%)) expanded to the second order 
in vortex position is: 

-tr(GoAff) = -tr[G (r v V v )ff ]~tr [G (r v V v ) 2 ff ] . 

In absence of disorder, translational symmetry requires 
that there be no linear terms in r v (r) in the vortex action. 
This is explicitly demonstrated in Appendix 1X1 At the 
second order in r v we have: 



itr[G (r v V v ) 2 ffo] 



(21) 



= - J drdV tr[G (r,T;r,T)(r v (r)V v ) 2 tf (r-)] . 

Since there is only one time variable in this integral, it 
produces a "trapping" potential ~ |t v (t)| 2 , which seems 
to trap the vortex at the origin even in absence of disor- 
der. This term will, therefore, be cancelled in the com- 
plete vortex action, as we will see below. 
We are left with the third term of l (T%|l : 



-tr(G Aff) 2 



(22) 



drd 2 ? 



dr'dVtr 



G (r,T;r',T')(r v (T')V v H (r')) 



G (r',r';r,r) (r v (r)V v ff (O) 



Let us utilize the absence of disorder and introduce a 
convenient notation for the following matrix elements: 



U n , n > = - J dW„(r) V v? Mr) 
= J d 2 r^ (r) VW(r) • 



(23) 



If we substitute l|19(l into the trace l|22(l . integrate out 
the time variables and partially sum over frequencies (de- 
tails of which are in Appendix [SJ , we obtain a simple 
expression: 



l -tv{G AH) 2 



(24) 



1 E 1 

9 i—i R 



f(e n )-f(e n ,) 



(e n - en') 2 | r v(w)t/ ni „'| 2 



where /(e) = 1/(1 + e^ e ) is the Fermi-Dirac distribu- 
tion function. Note that in the last expression u> already 
describes vortex motion - it is a "bosonic" Matsubara fre- 
quency, taking values 2kirT at finite temperatures. This 
expression is not zero at lu — 0, which implies its non- 
invariance under translations, even in absence of disorder. 
However, the trapping potential encountered earlier can 
now be exactly cancelled. A simple recipe, which is rig- 
orously derived in Appendix lAl is to subtract from i|24|) 
its zero-frequency part. This way we obtain the complete 
effective vortex action: 



S^tff ^ ft \\ iuj ^ n ~ £ "') 

(/( e «) ~ x 



e n i - iuj 



\r v (u)U nn , 



(25) 



The vortex action must be invariant under all lattice 
symmetries of the superconducting material. Transla- 
tional symmetry is already guaranteed, while the 90° de- 
gree rotation symmetry of the square lattice is imple- 
mented only after the contributions from all quasiparti- 
cle nodes are included. Symmetry under 90° rotations 
generates the same requirements on the vortex action at 
the second order in vortex displacement as the symmetry 
under any rotation. The allowed rotationally invariant 
terms are therefore |r v (u;)| 2 and iz(r*(ui) x (r v (u>)), and 
we may write the action as: 

5v = ^Eh H|rvH|2+F±H ^ (r:H x (rvH) 

UJ 

The first term describes vortex mass and dissipation, 
while the second term breaks the time-reversal symmetry 
and describes possible Magnus-type forces on the vortex 
arising from presence of quasiparticles. 



IV. QUASIPARTICLE HAMILTONIAN 

In this section we derive a convenient linearized form 
of the quasiparticle Hamiltonian that will be used to find 
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the effective vortex action. While our focus will be the 
Bogoliubov-de Gennes (BdG) description of quasiparti- 
cles, we note that similar results can also be obtained 
from the "staggered-flux" model of quasiparticles in d- 
wave superconductors. 

The BdG Hamiltonian for an extreme type-II d-wave 
superconductor is: 



-ffBdG = 



(p-f^) 2 _ E 



2m 



A* 



(p+tA) 



(27) 



where p is the momentum operator, A gauge field, Ef = 
h 2 kp/2m quasiparticle Fermi energy, and A gap operator 
given by: 



A 



J { Px , { Pv , A(r)}} - ^A(r){4, d y }$(r) 



h 2 k 2 F 

Mr) 
h 2 k 2 F 



Px 



hd x ${r) 



>Py 



hd y <P(r) 



(28) 



A(r) is the superconducting center-of-mass complex gap 
function, with phase $(r), and amplitude Ao assumed to 
be uniform (except at vortex cores). The braces denote 
anticommutators defined as {a, b} — ^(ab + ba). This 
Hamiltonian formally has d xy symmetry, which is related 
to d x 2_ y 2 by a rotation. 

We follow the approach from Ref. |3Jfl and 122. It was 
developed for a vortex lattice, but we will use it in the 
limit of infinite inter-vortex separation. A vortex lattice 
can be partitioned into two sublattices, A and B. Let 
$A and <1>b be superfluid phases that originate only from 
the vortices on A and B sublattices respectively, so that 
$(r) = <&A(r) + < &b(^)- Define: 



mv A {r) 
mvB(r) 

a(r) 
v(r) 



KV^ A (r) - -A(r)) 
c I 



W$ B W - -A(r) 

\ c 

mv A {r) — mvB(r) 
2 

v A (r) + v B (r) 



(29) 



In the extreme type-II limit the London penetration 
depth is large, and the magnetic field may appear prac- 
tically uniform across inter-vortex distances. In such cir- 
cumstances one can neglect contribution of the gauge 
field to the total flux of mv on A and B sublattices and 
keep only the 2tt contribution of the superfluid order pa- 
rameter. Then, the fields a(r) and v(r) look like vector 
potentials of a 7r-flux vortex lattice, where the flux al- 
ters sign between A and B sublattices in o(r), but not 
in v(r), which is the superfluid velocity. 

Now perform a unitary Franz- Tesanovic 
transformation 31 U on the Hamiltonian l|27(l : 



U 



-i<S> A (r) 








(30) 




D 

(p— mn B )' 
2m 



E, 



with the transformed gap operator: 



D = 



Ap 

h 2 k 2 F 



{p x + a x ,p y + a y } . 



(31) 



(32) 



Finally, this Hamiltonian can be linearized near any of 
the four nodal points: 



near ± hkpx 



(33) 



H = ±[vp(p x + a x )a z + v/±{p y + a y )a x + 7nvpv x I] 

near ± HkFy '■ 
H = ± [v F (p y + a y )a z + v A {p x + a x )a x + mv F v y I] , 

where a x ,a v ,a z are Pauli matrices, and / is the unit 
two- by- two matrix. The last term in these Hamiltonians, 
proportional to /, is the Doppler shift of quasiparticle 
energy due to the background superfluid flow. 

Before we take the limit of infinite inter- vortex separa- 
tion, it is important to understand how the magnetic field 
varies with distance from a vortex core. A more accurate 
expression for the superfluid velocity v(r) is obtained by 
solving the conventional London equation (Appendix of 
Ref. El): 



v(r) 



2nh 
m 



d k iky. z 
(2tt) 2 k 2 



1 



1 + X 2 k 2 



C k e lkr 



where A is the London penetration depth, and: 



-ikr- 



-ikrf 



Ck = 



(34) 



(35) 



Note that the field a(r) does not depend on A; all of 
A dependence is contained in the Doppler shift of 
The superfluid velocity v(r) decays as 1/r with distance 
from the vortex core when A — > oo, and faster for finite 
A. 



V. ANALYTIC RESULTS IN ABSENCE OF THE 
DOPPLER SHIFT 

The expression (I25f) for the effective vortex action, and 
the description of quasiparticles given by (|33|l provide us 
with all ingredients we need to study vortex dynamics. 

We will follow two routes to the study of this dynamics, 
and the reader may examine them in either order. 

A perturbative computation is presented in Section lVll 
and the reader can jump ahead to that section now. 
There we begin with an alternative formulation22iSi2ii2^ 
of the same theory using auxiliary Chern-Simons gauge 
fields. This allows a straightforward consideration of 



7 



the structure of the (uncontrolled) perturbation theory, 
which leads to the results in Eq. (JSJ. 

In the second route, we use a non-perturbative analy- 
sis, and show that certain physical effects can be evalu- 
ated exactly. This route proceeds in two stages. First, 
we will ignore the Doppler shift of quasiparticle energies 
in the present section, and this will allow us to analyti- 
cally derive the vortex action. Such analytic results will 
elucidate interesting physics, and serve as a framework 
for the second stage (section IVllfl in which we include 
the Doppler shift. Calculations of the Doppler shift ef- 
fects must be done numerically, but fortunately, no qual- 
itative modifications of the analytic partial results are 
found. All these results are also consistent with Eq. 10. 



A. Transition matrix elements 

The main problem is to diagonalize the linearized BdG 
Hamiltonian (|33[) and calculate the matrix elements 
U n}n ' given by (|23|l . which need to be substituted into 
the vortex action (|25|l . The matrix elements V ' n ,n' can 
be evaluated exactly in absence of the Doppler shift. For 
now we focus on the node k = hkpx, and postpone dis- 
cussion of the other nodes for later. We start by solving 
the Schrodingcr equation (|20|) for the Hamiltonian: 



H _ | v F (p x + a x ) VA(p v + a y ) 
v&.{py + a y ) -v F (p 



(36) 



We will approximate the "gauge" field a(r) to that of 
a 7r-fiux vortex completely localized at the origin. Even 
though the vortex core will appear infinitely small, the 
approximation is validated by the fact that there are no 
strictly bound quasiparticle states in a d-wave vortex, so 
that the vortex dynamics is determined only by extended 
states. 

We will determine the spectrum of H following the 
analysis of Ref. |33- As they showed, microscopic details 
of the vortex core can be captured in a single parame- 
ter. The Hamiltonian l|3"6T) could be easily diagonalized 



if it were isotropic. Therefore, we make it isotropic by 
applying a gauge transformation that deforms the vor- 
tex gauge field a(r) to a suitable "elliptical" shape, and 
then rescale the coordinates in such a way that v F and 
i>a become equal to unity and the "elliptical" gauge field 
is reshaped to the isotropic form: 



— yx + xy 
2(x 2 +y 2 ) 



(37) 



The coordinate rescaling goes as follows (the original co- 
ordinates are "primed"): 



w x = v F w x ,w y = V A W y 
d 2 r' = VFVA<i 2 r 

v' = #■ + y- #■ 

H = (p x + a x )a z + (p y + a y )a x 



vectors 

Jacobian 

gradient 

wavefunctions 

Hamiltonian 



It is convenient to apply a unitary transformation to the 
rescaled Hamiltonian and write it as: 



H 



{Px + a x ) + i(p y + a y ) 



(Px + a x ) - i(p y + a y ) 



(38) 

where a quasiparticle mass m has been added "by hand" . 
In a d-wave superconductor m = 0, and eventually we set 
m to zero, but formally a non-zero mass is needed to inte- 
grate out quasiparticles in the path-integral and perform 
summations over frequencies such as (|A3|) . Switching to 
the polar coordinates: 



Px = cos^-*|.)-^(-^ a x 
the Schrodingcr equation takes form: 



sin (p 

~ 2r~ 

COS (p 



(39) 



-i<f> 



% ^ 2 



m 

r 



u(r) 
v(r) 



= E 



u(r) 
v(r) 



By substituting u{r) = e l{ - 1 1 ^u(r) and v(r) = e ll ^v(r) the solutions are found to be: 

ip q ,i,k(r, 0) 



4ir\E\ 

VE-m ■ J_ l+ i (fcr)e^- 1 ) 
—iq\/E + m ■ J_ ; _ i (fcr) 



y/E-m- J ; _i(fcr)e i( '~ 1)0 
iqy/E + m- j ;+ i (fcr)e 1 ^ 



(40) 



(41) 



I < 



I > 
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States are characterized by the quantum numbers n = 
(q, I, k): "charge" q = ±1 (distinguishes particles-like and 
hole-like states), angular momentum I S Z and radial 
wavevector k > 0. Energy is E — qy/k 2 + rri 2 . Ji(kr) are 



Bessel functions of the first kind, 
are normalized as: 



These wavefunctions 



The zero angular momentum channel requires special at- 
tention. Square integrability of the wavefunctions re- 
quires only that: 



^q,o,h{r, (j>) 



k 



Att\E\ 



(42) 



yE — m ■ J_i(kr)e 

2 

iq\J E + m ■ Ji (kr) 



-i(f> 



VE-m- Ji (fcr)e-^ 

2 v ' 

—iq\jE + m ■ J_ i (kr) 



I = 



where is a parameter with arbitrary value within [0, 7r). 
This parameter is physically obtained from boundary 
conditions at the vortex core only when the core is not 
idealized by a point, but regarded as a finite region in 
space. Microscopic details of the vortex core enter vor- 
tex dynamics only through this parameter. In the fol- 
lowing we assume that all of the flux is contained in a 
small cylindrical shell: this yields 6 = for a vortex, and 
9 = -jj for an anti-vorte3*2^. Other choices do not intro- 
duce qualitative changes in the results that follow, but 
they may give rise to some quantitative corrections, to 
the vortex mass for example, which we briefly discuss in 
the Appendix El 

The presence of a vortex significantly affects quantum 
motion of quasiparticles, as is clearly seen from the lo- 
cal density of states (LDOS) expressed in the rescaled 
coordinates (for m = 0): 



At small distances from the vortex core LDOS diverges 
as 1/r due to the zero angular momentum state l|42p. 
and there are also certain oscillations. This behavior was 
studied earlier, and with addition of the Doppler shift it 
only acquires anisotropy, but does not qualitatively de- 
viate from 1/r dependence^. In realistic circumstances 
this divergence is cut-off by the finite core size. Note that 
the LDOS peak at r = does not correspond to bound 
states in purely e?-wave vortex cores (bound states could 
emerge only if a secondary order parameter, such as s- 
wave, were present in vicinity of the vortex). In fact, a 
careful analysis shows that the Hamiltonian (|36fl docs 
not support even a bound state exactly at E = 02£. Far 
away from the vortex core LDOS becomes constant and 
linear in energy as expected. 

We can now calculate the matrix elements V '„y in 
(|23|l . This calculation can be done directly in the rescaled 
coordinate system, with the states (|41fl and (|42|l that 
were obtained after a gauge transformation, rescaling and 
a unitary transformation: 



U, 



i) 



y d 



d^t iW (-- + -J-)^(, ) .(44) 



Here rii denotes (gi, h,ki). 

The calculation of (|44|1 is pretty tedious, so that we 
only sketch it in the Appendix El A closed form for 
arbitrary quasiparticle mass m is obtained by algebraic 
manipulations and elementary integrations, with help of 
Bessel function identities. It turns out that no singulari- 
ties are introduced in the end by taking the m — > limit, 
so that in the following we will consider only the case of 
interest m = 0. Since U nit n 2 are essentially transition 
matrix elements of the momentum operator, the only al- 
lowed transitions are between states whose angular mo- 
mentum is different by one. Defining a = I2 — h = ±1, 
we obtain for the node hkpx and zero quasiparticle mass: 



2ir 2 r 
1 

2ir 2 r 
ill 

2n 



cos(2|e|r) , £ g ^ 



7T 

1=0 

\e\r < 1 
\e\r > 1 



U r , 



(43) where 



:e /fe-?0 ( cr—+i— I U, 



x 

vf va 



{qi,l-a,k 1 ),(q 2 d,k 2 ) ) 

(45) 



U, 



(q 1 ,l-o-,ki),(q 2 ,l,k 2 ) 



al- 



-A^k 2 8{E 1 - E 2 ) - C a , quq2 



al- 



2<yqiq 2 E1+E3 
ix E\ — E 2 



1 r 



5S?e(«r(fe-fa)) 



r 



2+1 
2 



(46) 



Here, 5(x) is the Dirac delta- function, and: 
6(z) = 




a 



(12 



(7 
(7 
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Quasiparticle energies are now Ei = qiki. The expres- 
sion l|46l) defines infinities at Ei = E 2 completely, ex- 
cept when I = (a + l)/2. In that case, however, it can 
be shown by finite-size regularization that U at E\ = E 2 
diverges only logarithmically with the system size, which 
does not affect the important integrals over energies in 
the limit of infinite system size beyond what is given by 
the expression l|4t)|) . 



B. Summation over quasiparticle nodes and 
angular momenta 



The effective vortex action l|25|) depends on U ni ,n 2 
through: 



r v (uj)U ni ,n 2 \ 
\xr v (cu)\ 2 \ijr v (uj 



U, 



(qi,l-<?.k 1 ),(q 2 ,l,k 2 ) 



(47) 



VfVa 



-z(rl(w) x r v (w)) 



where we have used l|45|) . relevant for the quasiparticle 
node at p = hkpx. Note that this dependence on vortex 
coordinates is not isotropic; only final expressions with 
all quasiparticle nodes included become isotropic. It is 
easily checked from (|33f) that the linearized Hamiltonians 
at different nodes are related by unitary transformations 
and an exchange vf <-> «a as needed, in absence of the 
Doppler shift. Thus, in order to add contributions at all 
nodes we only need to symmetrize with respect to vf and 
t>A and multiply by the number of nodes: 



£ = l^.'-^.fa.^)l 2 x (48) 



Ha 
vfva 



Note that for our purposes we do not need to worry about 
inter- node scattering of quasiparticles. If we went back 
to the non-linear quasiparticle Hamiltonian and calcu- 
lated U ni ,n 2 for states that belong to different nodes, we 
could obtain singular behavior only when the momentum 
transfer is ~ Hkp. Such transitions are not stimulated by 
a slowly moving vortex, and cannot play a significant role 
in low-energy vortex dynamics. 

Owing to the degeneracy of the quasiparticle energies 
with respect to angular momentum, it is possible to an- 
alytically carry out summation over angular momentum 
channels in (|25|) . that is summation over I and a. Vor- 
tex mass and possible dissipation are determined by the 
following sum, which is just a geometric series that we 



evaluate from l|46[) : 

|t / (gi,i-cr,fc 1 ),( ?2 ,i,fc 2 )| 

Z,er 

E ± +E 2 



{Ex - E 2 y 



(E 1+ E 2 ) 2 (E 2 + E 2 ) 
\El-El\^{El,El) 

2 fh-k 2 
log 



1 / 1 



4U? + H )l0S 



k\ — k 2 
ki + k 2 



k\k 2 

2 



(49) 



Here, we have dropped the Dirac delta- functions in 146|l . 
since in l|25|l there is a factor that becomes zero when 
Ei = E 2 (safe to do only in absence of the Doppler shift). 
For the same reason the other possible divergences at 
k\ = k 2 , in (|46(l or in the sum over angular momenta turn 
out not be dangerous: they actually appear only when 
qi = q 2 as well, that is only when E\ = E 2 (particle- 
particle or hole-hole transitions). 

Similarly, we find the sum that determines "transversal 
dynamics" : 

yi " \U(gi,l-a,ki),(q2,lM)\ = 
l,tr 

mS^^M^I) 



Ei+E 2 



e 1 -e 2 \e 1 e 2 J B \k 1 + k 2 



fcl 



ki + k 2 



(50) 



C. Vortex mass at zero temperature 

This subsection collects our results above and com- 
putes the contribution of the 7r-flux to the vortex mass. 
We will find that this contribution is finite. 

The vortex mass is determined by the "longitudinal" 
part of the vortex action (|26|l . Thus we will restrict 
our attention to the sum (|49|l and substitute it into the 
vortex action Ij25(l . At zero temperature the difference 
between Fermi-Dirac functions f{£ ni )—f(^n 2 ) allows only 
transitions between states whose energies have opposite 
signs. This means that only particle-hole transitions will 
be allowed. Using the invariance of (|49|) under exchange 
Ei <-> E 2l and symmetrizing the "longitudinal" part of 
(|25|l with respect to its complex conjugate, we have: 

c f duj 1 at, Iaw u 2 (E 2 -Ei) 

5v|1 = J^J m J dE \E 2 -E 1 )* + ^ X 

-A 



(51) 



nodes l.a 



Energy cut-off is A. If slow vortex dynamics reduces to 
vortex inertia, then the leading low-frequency term in the 
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"longitudinal" vortex action is: 

S vn -- 



2tt 2 

where m v is the vortex mass: 



dwm v w 2 2 



(52) 



(53) 



-A 



_ rp l^(<?l, i-<r,fci),(q 2 ,i,' 



Note that this expression is obtained by simply setting 
u> = in the denominator of 151(1 . which is justified 
only if m v turns out to be finite. In absence of thermal 
fluctuations m v could in principle be infra-red divergent 
only due to transitions at E\ = E 2 = 0. However, only 
the bulk density of states, which shapes such transitions, 
can be responsible for infra-red divergences. In the bulk, 
far away from the vortex core, LDOS is linear in energy 
(see (@2J)) and hence suppresses transitions at E = 0. 
This is formally reflected in the sum l)49|) by absence of 
the singularity at E\ = + , E2 = 0~. Then, without help 
from the transition matrix elements, the extra factor of 
E2 — Ei in l(53|l is innocuous and fully integrable. Note 
that the 1/r behavior of LDOS close to the vortex core 
matters, but not enough to cause infra-red divergency of 
to v ; it produces the second term in 1)49(1 . including the 
logarithms, but no singularity at E = 0. 

Careful integration, presented in Appendix [U] shows 
that indeed the vortex mass is finite: 



D. Dissipation at finite temperatures 

At finite temperatures the "longitudinal" part of the 
vortex action is: 

A A 



-A -A 



uj 2 {E 2 - Ej) 
{E 2 -E 1 ) 2 +^ 



^^|r v (o;)f/ ni , n2 | 2 .(56) 



nodes /,i 



The difference between Fermi-Dirac functions f(Ei) — 
f(E 2 ) does not strictly prohibit transitions between two 
particle-like or two hole-like states at T > 0. Therefore, 
singularities in which the factors E 2 — E\ become zero 
at finite energies will be included in the vortex action. 
We can no longer naively neglect to 2 in the denominator, 
even for small frequencies. 

Consider the "longitudinal correlation" Fii(a;) defined 
by At small frequencies we expect that small values 

of E 2 — Ei will matter the most, so that we approximate: 



f(E 1 )-f(E 2 ) » (E1-E2) 



df_ 

dE 



\E 1 -E 2 \<€T. 



With a change of variables AE = E 2 - Ex , 2E = E 2 + E x 
we have at small frequencies: 



1/1 1 



T A 

dA£ / dE 



-T 



-A 



df(E)\ 
^~dE~) X 



lu AE 



' X! \ U (n,l-<rM),(q2,i,k2)\' 



(57) 



0.05 



4- + 4- 1 a 

Vp V 2 A 



(54) 



Therefore, vortex dynamics at zero temperature is en- 
tirely specified by inertia. The cut-off scale A is roughly 
given by the maximum superconducting gap Ao- We can 
write: 



0.05 a D 



1 

OLD 



(55) 



where m e is the electron mass, and old = vp/v&. In d- 
wave superconductors an ~ 1 4- 10 2 , so that the purely 
quantum effects contribute roughly an electron mass to 
the mass of a vortex. 

In various circumstances, which are not the focus of 
this paper, the bulk LDOS may not vanish at zero en- 
ergy. For example, this can be expected in presence of 
disorder, or Zeeman splitting in magnetic field. In such 
situations the vortex mass calculated from the expression 
(|53|l could become infra-red divergent, and it would be 
more appropriate to go back to the full action l]5ip. 



We do not care about \AE\ > T, because for sufficiently 
large and finite AE we can safely expand in powers of u> 2 
and obtain only the high-frequency corrections to vortex 
inertia. Our goal here is to elucidate sub-quadratic pow- 
ers of frequency. Similarly, it will be clear shortly that 
the precise boundaries of AE do not matter after analytic 
continuation. The sum over angular momenta is given by 
(|49|l . Its "most divergent" part (2 nd term), which comes 
from the particle-particle and hole-hole transitions that 
involve the zero angular momentum channel, dominates 
the small frequency behavior: 



1 



&AE / dE 



w 2 AE 2 
AE 2 + lj- 

7T 

6 



32 



-T 

E 2 



-A 



AE 2 



~dE~) X 

(58) 



After analytic continuation to real frequency, this has 
the "Ohmic" dissipation form of the Bardeen-Stephen 
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viscous drag. The ellipses denote corrections that become 
more noticeable as the frequency grows. The first such 
correction is also interesting and comes from particle- 
particle and hole-hole transitions that do not involve the 
zero angular momentum channel (1 st term of (|49|l % ): 



T A 

AW - + J dAE J dE 

-T -A 



df(E) \ 
dE ) 



oj 2 AE 2 



-A 
'2/0 J?2 



AE\2E Z 



2 



) 



AE 2 +lo 2 \2EAE\{E 2 



AE 2 
4 



\EAE\) 



log(2) 



A 



Tuj 2 log 1 



^ 2 



This could be interpreted as a frequency-dependent cor- 
rection to the vortex mass. With such an interpretation 
the vortex mass at rest would be logarithmically diver- 
gent. However, we feel that it is more natural to as- 
sociate only the strictly quadratic frequency behavior to 
the vortex mass, and regard everything that dominates it 
at small frequencies as various forms of dissipation. All 
other corrections are (D(ui 2 ), and hence contribute the 
normal vortex mass. 

In conclusion, dissipation occurs at finite tempera- 
tures, its strength being a power law in temperature due 
to the gapless quasiparticle nodes. 



E. Transversal dynamics 

The "transversal correlation" F±(lj) defined in l|26() 
describes aspects of vortex dynamics that are associated 
with the Magnus force. We argue below that the nodal 
quasiparticles do not give rise to transversal forces on a 
vortex at any temperature. 

The expression for "transversal correlations" is: 

A A 

F±(f») = / d£i J dE 2 (f(E 1 )-f(E 2 ))x 

-A -A 

-iu(E 2 -Ei) 



E 2 — E\ — iuj 



,v\U( qu i-<jM),(<i2,lM)\' 



(59) 



The sum over angular momenta is given by i)50[l: its 
crucial property as a function of E\ and E 2 is that it 
changes sign under the change of variables E\ — > — E 2 
and E 2 — > — E\. If we apply this change of variables 
to F±(uj), and add the result to the original expression 
above, we obtain: 

A A 

2F ± (u)oc J dE 1 J dE 2 J2v\U( qi ,i-a, kl ),(<i2,iM)\ 2 x 

-iuj(E 2 -Et) 



La 



{f{E l )+f{-E l )-f{E 2 )-f{-E 2 )) 



E 2 — Ei 



The final conclusion follows from the identity f(E) + 
f(-E) = 1. 

The physical reason for complete absence of 
quasiparticle-induced transversal forces is related 
to the fact that there are no circulating currents of 
quasiparticles in thermal equilibrium. This connection 
has been made at small frequencies in Ref. At small 
frequencies it is safe to neglect u> in the denominator 
of H59|) . and then cancel out the factors of E 2 — E\. 
What remains is an integral that can be represented 
as a pure trace of a quantum operator, and this trace 
can be evaluated in real space. Using the definition of 
transition matrix elements (|23[1 . one finds: 



F±(oj) oc — ioj J d 2 r£(V x j) = —iuj<j)jds 



(60) 



where j is the net current density of quasiparticles, de- 
termined by wavcfunctions (|41|) . 142|) . and temperature. 
The contour integral above is taken on a loop of infinite 
radius. In thermal equilibrium quasiparticles are at rest 
with respect to the substrate, so that j = 0. 



VI. PERTURB ATIVE COMPUTATION 

This section reviews an alternative 
formulation 2 ^ 2 ^ 2 ^ 2 ^ of the quantum theory of vor- 
tices, fermionic quasiparticles, and plasmons which 
appears in Sections [I] and IIVI In this approach, the 
complete action for the vortices, plasmons, and nodal 
quasiparticles is S v + 1S4 + + S cs specified in Eqs. 161|l . 
0, (|63() . and (|64|l . This theory is closely related to a 
continuum limit of the theory considered by Balents and 
Fisher m& 

The new feature of this approach is that the Berry 
phase terms are mediated by auxiliary Chern- Simons 
gauge fields. Such a formalism readily allows a pertur- 
bative analysis of the couplings between the elementary 
excitations using standard techniques. It must be kept 
in mind, however, that there is no small parameter which 
justifies such a perturbative analysis. 

We will see below that the perturbative analysis of the 
Berry phase terms lead to conclusions consistent with the 
exact evaluation already presented in Section Em- 
powered by this success, we will examine the influence 
of the Doppler shift term using perturbation theory in 
Section IVIAl 

We begin by extending the vortex action S v in Eq. (JTJ 
by coupling it in addition to a second U(l) gauge field 



S v = 



dr 



2 



dr v (r) 
dr 



+ i^.i(r T (r)) 



+ % 



J d 2 rdT{A f ,{r,r)+a ll (r,T)}J vll (r,T) 



(61) 
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This is to facilitate the Berry phase coupling to the 
fcrmionic quasiparticles, as we now describe. 

The fcrmionic quasiparticle Hamiltonian is discussed in 
SectionslTland lTVl and we work with a single node appear- 
ing in the first equation in Eq. (|33|l . Let this Hamiltonian 
act on the Nambu spinor We define W = — i^a y , in- 
troduce the Dirac gamma matrices 7 M = (— a v , <r x , er z ), 
rescale coordinates and gauge fields to obtain isotropic 
velocities as discussed in Section [Vj and thence obtain 
the continuum action for the fcrmions 



d 2 rdr 



(62) 

The last Doppler-shift term is a coupling to the superflow 
fluctuations, which we have represented in this paper by 
the U(l) gauge field A^. Extracting the explicit con- 
nection between the gradient of the phase of the order 
parameter and the 'electric' field associated with _4 M , we 
may rewrite this term to obtain the fermionic action in 
its final form 



d rdr 



IV p 



+ — ^tf 7 °* (OyAr - d T Ay 



(63) 



Finally, we need to tie the gauge field a M to the vor- 
tices, as was done explicitly in Eq. (|37|l . In the present 
formalism, this is done conveniently with a Chern-Simons 
term: 



S cs = - / d 2 rdT€ tl i,\a fJ ,d u ax 



(64) 



We have now written all the required terms coupling 
together the vortices, fermionic quasiparticles, and plas- 
mons in a d-wave superconductor. The complete action 
is S v + S^ + Sty + S cs specified in Eqs. (|£T)l . ©, (JHSJl, 
and 

At this point, it is useful to describe how the results of 
Sections IIVI IVl and IVIII relate to the present formalism. 
The Hamiltonian in Eq. (|33|1 is obtained by exactly inte- 
grating out the gauge fields a M , a^, and only the scalar 
potential A T in the Coulomb gauge V • A = 0. The 
retardation effects contained in A fluctuations are ne- 
glected. The computations of SectionlVland lvTTl desCTibe 
the consequences of subsequently exactly integrating out 
the fermions. 

The present section proceeds by integrating out de- 
grees of freedom in the opposite order. We first inte- 
grate out the fermionic fields \& exactly, and then study 
the perturbative consequences of subsequently integrat- 
ing out the a M , a M and A^. As we have already seen in 
Section [n] such an approach readily allows accounting 
for the A fluctuations. 

We will first integrate out the fermions VE", and study its 
influence on the and a M gauge fields, and then inves- 
tigate the resulting modifications to the vortex/plasmon 



action in Eqs. JTJ and J3J). We defer consideration of 
the Doppler shift term, and the resulting modification 
of plasmon (i.e. A^) fluctuations to the following Sec- 
tion eei 

Integrating out \& in Eq. l|63l) also yields an effective 
action Trln[7 M (<9 M — ia^)\ for a^. This involves terms to 
all orders in a M and is not easy to work with. To obtain a 
first understanding of its influence, we expand the func- 
tional determinant to second order in a^. Subsequently 
we integrate out the a M , while accounting for the Chern- 
Simons term in S cs . This yields an action for the gauge 



field a M of the form 



c f d^p f P P 

Sa= 2j 8^'^ ( S ^-^ 



(65) 



where c is a coupling constant. Now we can examine the 
effect of fluctuations on the vortices just as in Sec- 
tion [H] The action S a is similar to that in Eq. (JHJ but 
with an extra power of p. Consequently, we obtain results 
for the effective vortex action as in Eqs. I|llfl and l|12l) but 
with an extra power of 1/r in the integrand. The r inte- 
grals are now linearly convergent in the infrared, and we 
obtain then a finite renormalization of the vortex mass. 
This structure is entirely consistent with the exact results 
of Section [3 where we found vortex mass contributions 
expressed in terms of integrals which were also linearly 
convergent in the infrared. 

More explicitly, in frequency and momentum space we 
obtain an expression which is just as in Eq. @, but with 
the (u) 2 + k 2 ) in the denominator replaced by ^/uj 2 + k 2 . 
Expanding the resulting expression in powers of r v (u>) we 
obtain a result in the form of Eq. with a contribution 



1 fd 2 k {k 2 -kVWT^+2m 2 ) 

= S= / 7Z2 „ m , s (66) 



25 / 4tt 2 



2V k 2 +uj 2 



It is now easily seen that the result of the momentum 
integral (and after reverting to physical units by re- 
inserting factors of the velocities) has the form of Eq. JSJ, 
with F% a universal function of the velocity ratio. In other 
words, we obtain a finite non- universal renormalization 
of m v , and a universal "sub-Ohmic" dissipation. 



A. Doppler shift 

We are interested here in the influence of the last 
Doppler-shift coupling in Eq. (|63|l between the nodal 
fermions and the superflow (which is the dual 'electric 
field'). At first glance, it appears that this term is in- 
nocuous: integrating out the fermions leads to an ultra- 
violet finite fermion loop, and this contributes a term 
proportional to the square of the A^ 'electric' field with 
a finite co-efficient. Comparing with Eq. J3J we conclude 
that this is merely a renormalization of the superfluid 
stiffness p s which controls the plasmon fluctuations. 
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However, for completeness, we do have to also consider 
the subleading non-analytic momentum and frequency 
dependence of the fermion loop, which could lead to sin- 
gular corrections to the vortex action. We will see below 
that this does not happen, and that such terms lead only 
to a finite rcnormalization of the vortex mass. 

Let us write the electric field £ x — d T A x — d x A T , and 
similarly for £ y . Then the term 



to leading order in /C, we obtain the form in Eq. Q with 



in 2 p s 



d 2 r \ dr(£ 2 x +£ 2 ) 



in Eq. @ is replaced by 

1 f d 2 k f duj ' 
4^2 J 2^ 

fC(k y , k x ,oj) 



in 2 p s 



1 + 



Ps 



K,(k x , kyju) 

Ps 

\£ x (k,uj)\< 



(67) 

\£ y (k,cu)\ 2 
(68) 



We compute the universal correction K, at one loop order, 
using the original unsealed units in the Hamiltonian in 
Eq. (|33l) . This leads to (after including contributions 
from all nodes) 



IC(k x ,k y ,uj) 



d 2 p de 



Tr 



2 J Air 2 2tt 

x [-i(e + w) + v F (p x + k x )a z + v A (p y + k v )a x ]~ 

x [— ie + v F p x a z + vi±p y o x Y x 

2 f d 2 p de r . 2 . 

= v F ^-2— [-e{e + uj)+v F {p x + k x )p x 

+ v l(Py + k v)Py] V 2 + v FPl + v APy] 1 

x [(e + lu) 2 + v%(p x + k x ) 2 + v A (p y + ky) 2 ]' 1 



v f {lo 2 +3v 2 F k 2 x + iv\k 2 y ) 
64u A [uj 2 + v F kl + v 2 A kl\ 



,1/2 



(69) 



where dimensional regularization was used in evaluating 
the integrals (this is a rapid way of picking out the non- 
analytic piece, as the analytic contribution has already 
been absorbed by renormalizing p s ). 

Now we can proceed to a computation of the renormal- 
ization of the vortex action as in Section [H] except that 
the substitution in Eq. I)68|l has to be performed in the 
gauge field action in Eq. Q. As before, we work in the 
Coulomb gauge. We focus on the A T fluctuations, and 
the expression in Eq. I|13|) is now replaced by 



S v a = 2tt 2 p s I drdr' 



d 2 kdu 
K(k x ,k y ,uj)\ 2 



,— i(jj(r — t ) 



x e 



Ps 



ih- (t* v (r) — r v (r ')) 



K,(k y ,k x ,u) 



(70) 



d 2 k 1 



(2tt) 2 k 2 1 v 
kl{K.{k y , k x ,uj) - K(k y ,k x ,0))\ 



[kl(fC(k x , k y ,u) - fC(k x ,k y ,0)) 

(71) 



Expanding this result as usual in small r v (u>), and also 



It is interesting to note that all factors of p s have dropped 
out of the expression for -Fji(o;), and a universal expres- 
sion dependent only upon the velocities v F and t>A has 
been obtained. The momentum integral in Eq. I|71|) is 
also free of infrared divergences. However, an ultraviolet 
divergence is present, and then the result has the form 
claimed earlier in Eq. JSJ) with F2 a universal function of 
the velocity ratio. 

It can be verified in a similar manner that the fluctua- 
tions of the spatial component of the gauge field, A lead 
to a non-universal correction to the vortex mass. 



VII. INFLUENCE OF THE DOPPLER SHIFT: 
BEYOND PERTURBATION THEORY 

Section has provided insight into the vortex dynam- 
ics with neglected Dopplcr shift. Then we included the ef- 
fect of the Doppler shift in Section lYl Al but only in per- 
turbation theory. Now we go back to the full linearized 
quasiparticle Hamiltonian (|33J) and describe a compu- 
tation which includes the Doppler shift beyond pertur- 
bation theory. We will only focus on an isolated vortex 
at zero temperature (some interesting effects in presence 
of other vortices are discussed in the Appendix [DJ . In 
such circumstances, only the "longitudinal" dynamics is 
non-trivial. The main questions that we ask are how the 
vortex mass is modified, and whether dissipation or some 
anomalous dynamics emerges at small frequencies due to 
the Doppler shift. We will find that even after account- 
ing for the Doppler shift, the vortex mass remains finite, 
as anticipated in the perturbation theory in the previous 
section. 

It was hinted in Section[V]that infra-red divergent vor- 
tex mass, or various forms of dissipation, could emerge 
at zero temperature if the bulk LDOS were finite at zero 
energy. Generally, one would expect that the presence of 
a vortex does not alter the bulk quasiparticle LDOS. This 
was found to be true in absence of the Doppler shift, and 
should remain true when the Doppler shift is included, 
since it is a localized perturbation (at any finite Lon- 
don penetration depth). Then, the vanishing zero-energy 
bulk LDOS of an ideal d-wave superconductor leads to a 
finite quasiparticle contribution to the vortex mass. This 
scenario is intuitive, and needs to be verified by a rigorous 
calculation. We will also verify in this section whether 
the Doppler shift introduces dangerous changes in LDOS 
near the vortex core. The concern about possible anoma- 
lous vortex dynamics also comes form a different point of 
view. The anomalous dynamics was found in Section Ivl 
whenever particle-particle or hole-hole transitions were 
allowed. With the Doppler shift viewed as a perturba- 
tion, the exact eigenstates will be superpositions of states 
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(|41|l and (|42|l . so that the exact particle-like states will 
have "hole tails". Even the particle-hole transitions be- 
tween exact states, allowed at zero temperature, will in- 
corporate hole-hole transitions between the unperturbed 
states. 

In the following we will exactly diagonalize the Hamil- 
tonian on a sample with finite size R, and at the 

node p = hkpx; the other nodes are easily included on 
symmetry grounds. Only the worst case scenario of in- 
finite London penetration depth will be considered. We 
will use the states (|41|l and (|42|) as the basis states for 
our representation. The main reason for choosing these 
basis states is that the transition matrix elements be- 
tween them are analytically known. This way all singu- 
larities, which are difficult to handle numerically, will be 
captured inside analytical expressions. Therefore, even a 
rough information about wavefunctions in the thermody- 
namic limit R — > oo, obtained from numerical finite size 
scaling, might suffice to establish properties of the vortex 
action. 

At zero temperature we adapt the expression (|25() and 
write the longitudinal part of the vortex action: 

2A 

S vll = f^fdAe 



2tt J Ae 2 +uj 2 
o 



K(Ae)\r v (u)\* (72) 



TZ(Ae) cx J2 &{enM-e n ')S(e n - e n , - Ae)\U n , n , 



We will use labels n for the exact eigenstates, and e n for 
exact energies. The basis states will be labeled by their 
angular momentum I and energy E, with q = sign(_E) 
and k = \E\ understood. The exact transition matrix 
elements [J„y, defined by (|23[) . take the following form 
in the chosen representation: 



E 



J &E X J dE 2 (E 1 ,l 1 \n 1 )(n 2 \E 2 ,l 2 ) x 



(73) 



where the matrix elements U (E 1 ,i 1 ),(E 2 ,i2) between the 
basis states are given by (|4*5|) . 



A. Exact quasiparticle spectrum and states 

We carried out numerical diagonalizations with typi- 
cally Nk = 50-^100 discrete values for the radial wavevec- 
tor k. This merely sets the energy cut-off to A ~ NkR -1 , 
and provides a plenty of states with kR ^ 1 that we are 
interested in. The number of angular momentum chan- 
nels I used in calculations was varied between Ni = 9-^29. 
Ideally, Ni ~ AR is needed in order to include all possi- 
ble states with \E\ < A in the Hamiltonian representa- 
tion, but this would be too formidable. Fortunately, even 
though we worked with small Ni , the final results showed 



virtually no dependence on Ni , meaning that all interest- 
ing physics at low temperatures is set by small angular 
momenta. The finite size scaling can be done in two ways, 
by changing either R, or A. This is because there is no 
specific scale in the linearized Hamiltonian and thermo- 
dynamic limit (all coordinate dependence comes from the 
products kr) . Details of the Hamiltonian representation 
are shown in the Appendix IF1 

The numerically obtained spectrum appears to be gap- 
less. It also appears symmetric under "particle-hole" ex- 
change: for every state with energy e there is exactly one 
state with energy —e. This property is, actually, pro- 
tected by a symmetry of the full Hamiltonian. There is a 
unitary transformation that transforms the Hamiltonian 
H in (|33fl at any quasiparticle node to —H, and this 
transformation is simply a coordinate system rotation by 
180° degrees (<f) — > 4> + t m <|33[) after substitutions like 
(J2U). Therefore, the Doppler shift does not cause "spec- 
tral flow" , that is a shift of some hole-like energy levels 
(e < 0) into particle-like regions (e > 0) or the other 
way round. The situation would have been different if 
we considered a system of vortices, say a vortex lattice, 
that is not inversion symmetric. If the 180° rotation sym- 
metry were violated, "spectral flow" could happen. As 
analyzed in Appendix [DJ this could lead to dissipation 
and anomalous dynamics at zero temperature. 

An example of numerically obtained eigenstates is plot- 
ted in Figures EJa) and H£b). These figures show the 
amplitude of *f! n (E,l) — (E,l\n) for a randomly chosen 
eigenstate \n) of the full Hamiltonian as a function of E 
and I (which characterize the basis states - eigenstates 
in absence of the Doppler shift). In order to gain in- 
sight about what happens in the thermodynamic limit, 
the wavefunctions ^ n (E,l) were normalized to the sys- 
tem radius — instead of usual unity. The Figure Eta) 
compares two different system sizes R. It can be no- 
ticed that roughly at e ~ E there is a single point where 
the amplitude is l^l 2 cx R. This indicates that a Dirac 
peak may be developing in the thermodynamic limit at 
e = E. At all other points, the amplitudes \^>\ 2 coin- 
cide for the different system sizes, and there are "tails" 
of non-zero amplitude that spread away from e ~ E, de- 
scribable by finite analytic functions. Possible survival 
of the Dirac peak can be understood in a perturbative 
point of view where the strength of the Doppler shift 
Hjj is parametrized and the full Hamiltonian written as 
H = Hq + XHd ■ For A = 1 we would obtain the Hamil- 
tonian il-vil) that is being numerically diagonalized. If A 
were zero, then the wavefunction ^ would be only the 
Dirac peak by definition: ^/ e .i(E, I) — 6ij6(e — E), while 
for A > the amplitude would start leaking to other 
energies e ^ E, as was found numerically. From the scat- 
tering theory we would expect that even for A = 1 the 
Dirac peak survives at e = E; the Doppler shift is a suf- 
ficiently localized perturbation to the Hamiltonian that 
still has a singular gauge field. 
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(a) 




(b) 



FIG. 2: (a) Typical wavefunction amplitude as a function 
of energy and angular momentum for two system sizes, R ~ 
63 and R w 125. (b) Typical wavefunction amplitude as a 
function of energy and angular momentum for three numbers 
of angular momentum channels, iV; = 9, Ni = 13, and Ni 
1 7. Normalization is — , and old ~ 0.4. 
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FIG. 3: Function IZ(Ae) for ao = 10, calculated with differ- 
ent numbers Z max G {4, 6, 8} of included angular momentum 
channels: I G { — i max , . . . , imax}. Both the horizontal and ver- 
tical axes are expressed in units that scale as 1/R with system 
radius R. 



B. Vortex dynamics at small frequencies 

In order to obtain Ohmic dissipation, the function 
IZ(Ae) in (|72|l should depend on Ae as Ae -1 . Actually, 
Ae" with any v < is a candidate for anomalous be- 
havior that dominates over finite vortex inertia at small 
frequencies. The plot of IZ(Ae) evaluated numerically is 
shown in the Figure [3J There arc no indications that any 
anomalous dynamics could occur. The situation resem- 
bles very much the one that was encountered in absence 
of the Doppler shift at zero temperature. The function 
IZ(Ae) appears to be linear, so that its slope s from the 
Figure |3 determines the vortex mass: 



s(ap) 
8 



1 



(74) 



For «d = 10 this roughly equals 3.3 electron masses, 
which is a factor of 6.59 enhancement of the vortex mass 
to that obtained in absence of the Doppler shift. Clearly 
the Doppler shift alone provides most of the vortex mass, 
but otherwise seems not to alter the qualitative physics 
that was uncovered in Section^ In Figure 01we compare 
the functions IZ(Ae) for several values of the anisotropy 
factor qdi and plot the slope s(au)- The slope s(ao) 
does not change significantly above ao > 10, so that in 
this range the vortex mass is roughly proportional to the 
anisotropy factor cvrj. In the limit au — * the matrix 
elements of the Doppler shift vanish (see Appendix [FJ , so 
that the slope s(ar>) slowly converges to s w 8*0.05 = 0.4 
given by (J5SJ. 

Another indication that there should be no anomalous 
dynamics is found in the quasiparticle density of states. 
The total density of states at low energies in Figure [S] is 
easily obtained from the exact spectra, and it is indica- 
tive of the bulk LDOS. The bulk LDOS remains a linear- 
function of energy, and roughly equal to LDOS in absence 
of the Doppler shift, as expected. Close to the vortex 
core LDOS acquires anisotropy in the rescaled coordi- 
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FIG. 4: (a) Function TZ(Ae) for a D G 

{0.33,1,5,10,100,1000} (slope grows with Qd), calcu- 
lated with Z max = 4. Both the horizontal and vertical axes 
are expressed in units that scale as 1/R with system radius 
R. (b) Slope of IZ(Ae) (linear fit) as a function of the 
anisotropy ratio an- 



FIG. 5: (a) Total number of quasiparticle states in presence 
of a vortex. Also shown is a quadratic fit to the data, (b) 
Total density of states of the nodal quasiparticles in presence 
of a vortex. Dashed line shows the density of states in absence 
of the Doppler shift. Both plots have parameters «d = 10, 
'max = 14, and both the horizontal and vertical axes are ex- 
pressed in units that scale as 1/R with system radius R. 



nates, but does not deviate from the 1/r behavior. An 
asymptotic solution^ for the wavefunctions that domi- 
nate LDOS at r — > in presence of the Doppler shift (at 
the node hkpx) is: 

(/ cos (h \ ( cos ch \ \ 

cxcosl -yM -c 2 sm(^J \ 

where c\ and c-i are constants fixed by matching to the 
regions of large r. The transition matrix elements ill'-il) 
that involve these states and other normalizable states 
are always finite; the only possible concern are the "di- 
agonal" matrix elements, but they must be finite as well 
because the states at finite energies carry finite momenta. 

The last argument against dissipation and anomalous 
dynamics is presented in the Appendix [HJ There we 
study in detail analytic properties of the exact transi- 



tion matrix elements U nit n 2 . Making use of the known 
transition matrix elements H45[) between the basis states 
in absence of the Doppler shift, we ask the following gen- 
eral question: what should the exact eigenstates look like 
in this basis in order for dissipation to occur. The gen- 
eral conclusion is that dissipation would require from the 
exact Hamiltonian to either systematically shift or shuf- 
fle the energy levels of the basis states in vicinity of the 
zero energy. The first scenario is equivalent to an effec- 
tive chemical potential for the linearized quasiparticles, 
and it was discussed in the Appendix |Dj), while the 
second scenario could happen due to Zeeman splitting or 
disorder. In both cases the bulk LDOS would be finite. 
Neither of these scenarios seem to be realized in our situa- 
tion of interest, when we consider the quantum-coherent 
Doppler shift due to the circulating supercurrents of a 
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single vortex. 



VIII. DISCUSSION AND CONCLUSIONS 

We studied various effects introduced by nodal quasi- 
particles in clean d-wave superconductors on dynamics 
of isolated vortices in thermal equilibrium. These effects 
generally give dominant contributions to vortex mass 
and dissipation in fermionic superfluids, over the corre- 
sponding contributions of the superfluid order parame- 
ter. In d-wave superconductors the quasiparticles have 
gapless nodes, and this makes quantum-mechanical phe- 
nomena fundamentally important for vortex dynamics 
at low temperatures. Therefore, we used methods that 
carefully took into account the phase change of quasi- 
particle wavefunction upon encircling a vortex, resulting 
with the proper determination of quasiparticle density 
of states. Classical effects of the superfluid flow in the 
vicinity of a vortex were taken into account through the 
Doppler shift of quasiparticle Hamiltonian, and it was 
shown that they do not lead to qualitative changes, but 
give rise to important quantitative corrections to the vor- 
tex mass. Vortex dynamics was studied both perturba- 
tively and non-perturbatively, and the two approaches 
yielded the same results. The perturbation theory was 
formulated in a dual description of d-wave superconduc- 
tors, where the Cooper-pair fluctuations are handled as 
a U(l) gauge field that mediates interactions between 
"charged" particle-like vortices. The non-perturbative 
approach was focused on the small-amplitude response of 
a vortex to a weak oscillating force. Our findings are also 
in agreement with simple scaling arguments that regard 
the d-wave quasiparticles as a quantum critical system. 

Our main conclusions are the following: (1) at zero 
temperature the vortex mass is finite and completely 
specifies vortex dynamics, apart from a universal "sub- 
Ohmic" damping; (2) at finite temperatures vortex mo- 
tion becomes dissipated in an Ohmic fashion, with cer- 
tain logarithmic corrections. Furthermore, quasiparti- 
cles mediate new interactions between vortices (assumed 
to be far enough apart so that quasiparticle motion be- 
tween vortices is incoherent), leading to both longitudi- 
nal and transversal forces in addition to the usual "hy- 
drodynamic" forces. All parameters that characterize 
vortex dynamics depend quantitatively (and not quali- 
tatively) on the boundary conditions at the vortex core. 
We reached these conclusions under several assumptions. 
First, the vortex core was idealized by a point in the 
continuum limit. This was primarily motivated by the 
fact that in good d-wave superconductors the coherence 
length £ is not much larger than the Fermi wavelength 
kp 1 , so that the core is indeed small. Furthermore, the 
superconducting gap vanishes in the nodal momentum di- 
rections, leaving a route for quasiparticles inside the core 
to leak out. Both of these facts conspire against bound 
core states, which is the situation that our idealization 
of the core models well, and various other studies sup- 



port. The vortex mass is, then, determined by scattering 
of free quasiparticles from a moving vortex. Smallncss 
of the core renders quantum phenomena very important. 
At finite temperatures we ignored thermal fluctuations of 
the superfluid order parameter, but kept the quasiparti- 
cles in thermal equilibrium. Presumably, as long as tem- 
peratures are low and quasiparticle contributions to vor- 
tex dynamics dominant, such an approximation is valid. 
Effects of thermal non-equilibrium have not been taken 
into account, but they are expected not to be important 
in the circumstances that we consider (slow oscillatory 
vortex motion with very small amplitude) . 

As we noted in Section [I] there is a possible dis- 
agreement between our results and the semiclassical 
approachii*i2iiii4. The latter works find that the nodal 
quasiparticles induce an effective mass which diverges lin- 
early with the separation between vortices. In particu- 
lar, it does appear that direct application of our theory 
and the semiclassical theory leads to distinct predictions 
for a certain class of experiments. It has recently been 
argue d 16 ! 18 that the zero point quantum motion of vor- 
tices leads to modulations in the local density of states, 
and this was proposed as the explanation of recent STM 
experiments on the cuprate superconductors 38-39 . The 
spatial extent of the zero point motion, and hence den- 
sity modulations, is determined^ by the value of m v . (In- 
deed, this fact was the initial motivation for our study.) 
Consequently, m v can be determined, in principle, from 
STM measurements of the modulations. Our theory pre- 
dicts that such an m v should be independent of the field 
H, at small H, while the semiclassical theory predicts a 
divergence in m v as ~ 1/V7J. 

In the future, it is also possible that dynamic (THz 
scale) measurements of vortex dynamics will lead to 
more direct quantitative determinations of m v , and these 
should allow tests of the theoretical proposals. 

At the end, we make some comments about the inco- 
herent quasiparticle mediated interactions between vor- 
tices, explored in the Appendix[Dj They were modeled by 
an effective chemical potential /i in the linearized quasi- 
particle Hamiltonian, which describes the Doppler shift 
due to the supercurrents from a second nearby vortex. 
The interactions that we discussed in Appendix [D] are 
specific to d-wave superconductors (for similar effects to 
be observed in s-wave superconductors, /i would have to 
be larger than the gap, that is the two vortices practi- 
cally overlapping). The effects of these interactions would 
be visible in vortex-vortex scattering events: the vortices 
would lose some of their initial kinetic energy during scat- 
tering, and excite some low-energy extended quasiparti- 
cle states. Anisotropy of d-wave superconductors would 
be reflected in these interactions, so that the differential 
scattering cross-section, and the amount of kinetic energy 
lost, would depend on the initial direction of approach to- 
ward the point of collision. We also speculate that these 
interactions might play an important role in shaping the 
properties of vortex lattices (especially if pinning disorder 
ruins the vortex lattice inversion symmetry) . They could 
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lead to anisotropic stiffness and pinning of the vortex lat- 
tice orientation to that of the substrate. The discovered 
dissipation-looking response of a single vortex in pres- 
ence of another vortex could transform into significant 
and magnetic field dependent correction to the vortex 
mass in the equilibrium conditions of the vortex lattice. 
We leave these interesting issues for future work. 
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APPENDIX A: DETAILS OF THE VORTEX 
ACTION DERIVATION 

In this appendix we provide support for the section UTTI 
by explicitly evaluating traces from the expansion (|18f) . 
At the first order in vortex displacement we have: 

- tr[G (r v V v )#o] = 

= -J drd 2 r tr [G (r,r;r,T)(r v (r)V v )H (r)] 



drd 2 r ^ 



-y 



(Al) 



tr [tW^(r)r v (r)Vffo(r)] 

1 e iuj0+ 
drd r > — > (e, 

n uj 

te[(r v (r)VV„(r))4(r)] = °> 



where we have used lfH?|) . (|2U|) . and integration by parts. 
For integration purposes, in absence of disorder, we may 
replace the vortex position gradient V v with the negative 
normal gradient: V v — > —V. This is just a shift of the 
coordinate system that puts a displaced vortex back to 
the origin, and simplifies derivations. 

At the second order in vortex displacement we have 
two contributions. One comes from the second term in 



-itr [G (r v V v ) 2 #o] 
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„iujtT 
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In order to evaluate it, we sum over the frequencies: 
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Then, after integration by parts, some algebraic manipu- 
lation and making use of the fact that the wavefunctions 
ipn(r) form a complete set of states, the expression (1A2|) 
reduces to: 

\ (/( £ «) ^ /( £ «')) " en) J dr|r v (r)t/ n ,„f 

n, n' 

(A4) 

This has the form of a trapping potential. 

The other contribution to the second order in vortex 
displacement comes from the third term in l|18(l : 



itr(GoA^) 2 = 
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G {r',r';r,r) (r v (r)V v H (r)) 
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The matrix elements U 7un i are defined by 123|) . If we 
integrate out r and t' and use U n ',n — ~U nn i we obtain: 



itr (G „Ai/) 2 = \y^y 



-i(uj—uj')(t—t ) 



' (3 2 ^ {-iuj + e n )(-zu/ + e n t) 



\2|„ I. . , ,I\TT |2 



(e n -e n ') |r v (w - w')Z7 nj „/| J 



(A6) 



Let Aw = w — lo' . We can carry out summation over w 
using the formula: 



-y 



1 



P ^ (-iuj + ei)(~iuj + e 2 ) 
The trace (IA6I) becomes: 



f(ei) ~ f(e 2 ) 
ei - e 2 



(A7) 



1 



tr{G AH) 2 



-y-y 



f(e n )~ f(e n ,+iAu) 



n.n A 



€71 €-11 
2|„ ( A , ,\TT |2 



iAlu 



(e n -e n /)^|r v (Aw)Z7 n ,„/|" 



(A8) 



Note that Aw takes values 2irT x integer, so that /(e„< + 
iAw) is simply /(e„/). 

The vortex action at the second order in vortex dis- 
placement is obtained by adding l|A4(l and (|A8I) . In 
expression (|A4fl we perform Fourier transformation of 
r v (r) and integrate over r, while in expression l|A8(l we 
relabel Aw into w. The result is given by (|25|l . 
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APPENDIX B: NOTES ON THE DERIVATION 
OF t/ nil „ 2 IN ABSENCE OF THE DOPPLER 
SHIFT 

Here wc outline some steps that lead to the result (|43)) . 
We calculate (|44|l in the rescaled coordinate system upon 
substitution of l|41|l and i|42|) . together with the rep- 
resentation of gradient in the polar coordinates (|39|l . 
In order to treat all cases of angular momentum on the 
same footing in (|41|l and 142(1 , we label the indices of the 
Bessel functions in the upper and lower spinor component 
by V and I" respectively. Later we worry about depen- 
dence of V and I" on I. First we integrate out the polar 
angle <j>. This prohibits all transitions between states 
whose angular momenta do not differ by one. Hence, let 
a = I2 — h = il- Then, we use the following identities 
for the Bessel functions: 



9 = we can write: 



I' = 



-i + h 



1" 



i + h 



1 < 
1 > 

1 < 
1 > 



(B5) 



If this is substituted in i|B4|) . than it can be easily seen 
that the angular momentum indices of 7 never differ by 
more than two. This allows us to evaluate the 7 terms. 
If the angular momentum indices are equal, then: 



li,i{k 1 ,k 2 ) = —S(k 2 - fci) 

«2 



(B6) 



If the angular momentum indices differ by two, then use 
the first identity from IjBlfl to obtain: 



/ 

kr 
ld_ 

k d r 



Ji{kr) 
Ji{kr) 



(Ji-i(fcr) + Jj+i(fcr)) (Bl) 
(Ji-i(fcr) - Ji+i{kr)) 



li,i±2{ki,k 2 ) = 2 ^. ± — / drJi(k 1 r)Ji ±1 (k 2 r) 



li,i{kiik 2 ) 



(B7) 



in order to carry out differentiation over r and eliminate 
the 1 jr factors that originate from the azimuthal gradient 
components. We introduce notation: 



r yi 1 ,l 2 (h,k2) = J drrJ h (kir)Ji 2 (k 2 r) 




and 




y/Ei - m ) ^/E2~ 



A" = <?i<7 2 sign ( h - -j sign (l 2 - - I \ 




After some algebraic manipulations, we arrive at: 



tt , x y\k 2 

U ni ,n 2 = hid — X 

VF va I 4 



A 1 \\ 

A 1 

A" ( 1 



h-l 



(B2) 



(B3) 



(B4) 



H' 1 ,i' 2 -i(h,k2) + 



A" -1 



The integral above is a special case of Weber-Schafheitlin 
formulae, which yields simple expressions: 



1 



5(Jti 



(B8) 



k\ 



1+2 



7M +2 (fci,fc 2 ) - ^J^f|) B(A,- /,,: 



1 



5(fc 2 - fci) , 



where 0(x) is a step-function: 
0(aO = 




(B9) 



If Zi = or l 2 = then in some cases the two an- 
gular momentum indices of the 7 terms may differ by 
one. Then, however, exact expressions for the appropri- 
ate Bessel functions are simple enough to directly inte- 
grate out 7 and add their contributions to [7 ni! „ 2 . Fur- 
ther manipulations are still quite complex, and we do not 
include them here. We only note that it is safe at this 
stage to set m = (no singularities are introduced). The 
final result is l|45() . 



APPENDIX C: CALCULATION OF THE 
VORTEX MASS AT ZERO TEMPERATURE IN 
ABSENCE OF THE DOPPLER SHIFT 



for a = l 2 — l\ = ±1 (otherwise, U ni n2 = 0). Next, we 
evaluate 7; 1 ,i 2 (fci, £2)- Note that for the case of interest 



In this appendix we substitute (|49|) into l|53|l term by 
term, and explicitly calculate the integral for the vortex 
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mass in absence of the Doppler shift. The goal is to show 
that there are no infra-red divergences, and to provide 
an estimate of the vortex mass in terms of the electron 
mass. 

The first term of 149( 1 contributes the following to the 
vortex mass: 



/ J £2-ei |e 2 -e^|max(e 2 ,ef,) 
1 (fe - hf(kf + k 2 2 ) 



-A 
A A 



= / dfci J dk 


A fci 



' k 2 + k\ \k\ — kf |max(fc 2 , fc|) 



2 / dki / dfc : 



(h - k 2 )(k 2 + k 2 ) 

! (fcl + fc 2 ) 2 fc 2 

o o 
3(3-41og2)A = 0.682234A, 



where the constant of proportionality is: 



C 



(CI) 



easy to show that there are no infra-red divergences: 



Mi 3 ) 



A 



i 

f (x 1 -x 2 ) 2 
dxi / dx 2 - x 



7 b + 3 lo s : 





2 / Si 



Xi + x 2 
2 

X2 ' 



2A 

~~2 



tt/4 



Xi 4 

1 / cos c 



.j; 2 



log 



7T 



rdr 



£lX 2 

— sin 2c 



xi — x 2 
Xi + x 2 



sin 2(f) 



log 



2A 



o o 
1 — sin 2(f> 
1 + sin 2(f) 

tt/4 



r cos < 
1 



7T 



1 



sin 2 2(f) 
1 — sin 2(f) 



log 2 



sin c 
1 



sin 20 



cos cos ( 



sin < 



sin 2(f> 



log 



1 



1 — sin 2(f) 

1 + sin 2(f)) ^ sin 2 2(f) 
= -0.196462A , 



log 2 



1 + sin 2(f) 



1 — sin 2(f> 
1 + sin 2(f> 



where xi = fci/A = rcos</>, x 2 = fc 2 /A = rsinc6, and 
symmetry with respect to xi «-> x 2 was used. It is easy 
to check that the function of <j> inside the last integral is 
finite for every (f> E (0, 7r/4). The total vortex mass is: 

m v = + + m<® w 0.05 ( 4" + 4" ) A • ( C2 ) 



This contribution to the vortex mass includes effects of 
all particle-hole transitions between states that do not 
involve the zero angular momentum channel. No infra- 
red divergency has been encountered. The second term 
of Ij49(l includes partially the effects of transitions that 
involve the zero angular momentum channel. This term 
is the most divergent as E 2 — E% — > 0, and thus the most 
dangerous. However, it only gives a finite contribution 
to the vortex mass: 



oc /d £l /d e2 ^Af£l±£i) 

J J £2 - ei 7T 2 \ £ 1 _ £ 2/ 

-A 

A A 






A A-Afc/2 

Afe/2 

= ^ ^log2 - ij A = 0.313118A . 



The rest of (|49|) is hard to integrate exactly, but it is 



APPENDIX D: QUASIPARTICLE MEDIATED 
INTERACTIONS BETWEEN VORTICES 

In this appendix we will examine how the gapless 
quasiparticles of a d-wave superconductor mediate cer- 
tain interactions between vortices. Consider placing a 
second vortex near the vortex whose dynamics we want 
to study. We assume here that the second vortex is 
far enough so that the fermionic quasiparticles near the 
two vortices can be treated independently. This requires 
some external source of dissipation, which will render 
quasiparticle motion between the vortices incoherent. We 
assume in the remainder of this appendix that such a 
source of dissipation is present, and so our conclusions 
here are limited to this case. 

Recall the full linearized Hamiltonian for quasiparti- 
cles lj3l?|l . Presence of a second vortex affects the quasi- 
particles that move around the first vortex in two ways. 
First, the "gauge" field a(r) is changed, but this is not 
crucial because we can always "gauge" it away. Second, 
the superfluid velocity v(r) is changed, which modifies 
the Doppler shift term. If the distance between two vor- 
tices is large, the supercurrent due to the second vortex 
appears practically uniform in vicinity of the first vor- 
tex. The corresponding Doppler shift formally looks like 
a chemical potential y, presented to the linearized quasi- 
particles. This effective chemical potential depends on 
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nodal 
direction 



FIG. 6: One vortex in vicinity of another. Supercurrents due 
to the second vortex appear roughly uniform in the region 
occupied by the first vortex, resulting with a uniform Doppler 
shift of quasiparticle energies. 



the distance between vortices, and their orientation with 
respect to the nodal directions (/z oc v x , where v x is the 
component of superfluid velocity in the nodal direction 
x). 

In order to explore qualitative aspects of these interac- 
tions, we ignore complications arising from the vortex 
sclf-Doppler shift and non-uniformity of the supercur- 
rents. We simply add a chemical potential fil to the 
Hamiltonian ((36(1 . Effect of this is that all energies are 
shifted by \i in the vortex action. We restrict our atten- 
tion only to zero temperature, and assume /i > without 
loss of generality. Consider first the "longitudinal corre- 
lations" : 



\1 IV 

oc J &E X J dE 2 



m 2 (E 2 -E 1 ) 
(E 2 -E 1 ) 2 +lu 2 



Clearly, the anomalous dynamics at small frequencies and 
non-zero effective chemical potential is of the same kind 
as that found in the case of finite temperature. Only the 
energy scale that controls it has been changed. 

The remaining corrections are 0{u> 2 ), meaning that 
vortex inertia is also affected by the presence of another 
vortex. Physical consequences of all these effects could 
be directly observed in vortex-vortex scattering events 
as vortex kinetic energy loss to the quasiparticles. Also, 
dynamics of a vortex lattice should be affected. Due to 
quantum fluctuations of vortices the quasiparticle medi- 
ated forces would tend to pin the orientation of the vortex 
lattice to the substrate, and they would influence its stiff- 
ness. These effects would be even greater in absence of 
perfect inversion symmetry of the vortex lattice, caused 
for example by vortex pinning disorder. 

Now consider the "transversal correlations" : 

M A 

F ± (w) oc ( &Ej dE ^^j- El) x 



(D4) 



E 2 — E\ — iuj 

-A p, 

^2<7\U( qi J-aM),(q2,l,k 2 )\' 



La 



As shown in subsection IV El the symmetry under E\ —> 
—E 2 and E 2 — + — E\ would set F±(u) = 0. However, this 
symmetry is removed by finite \i. By substituting 1(50(1 
we find: 

F±(u) oc -ioj (/iA + 0(> 2 sign(X))) + C(w 3 ) . (D5) 

This result can be qualitatively understood by noting 
that the formal particle-hole symmetry of the linearized 
quasiparticle Hamiltonian is broken when /i 7^ 0. There 
is a net quasiparticle current density, formed by occupied 
states within the energy range < E < fi. This is possi- 
ble because the exact wavefunctions (141(1 and 1)42(1 are 
mixtures of a particle and hole that move with different 
angular velocities. 



(Dl) 



nodes /.- 



Similar to the case of finite temperatures and /1 = 0, cer- 
tain "dangerous" particle-particle transitions are allowed 
when [i > 0. They occur at E\ = E2 = fJ>, where the bulk 
LDOS is not zero. Following the procedure from subsec- 
tion [^^J we find that the transitions involving the zero 
angular momentum channel dominate at small frequen- 
cies: 

F\\(w) oc ^ 2 M , |«|«/*. (D2) 

The first correction at small frequencies comes from the 
transitions that do not involve the zero angular momen- 
tum channel: 

M oc fno 2 log (t^t) , H«/i- (D3) 
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APPENDIX E: OTHER VALUES OF THE CORE 
PARAMETER 6 

All calculations in Section W\ were restricted to bound- 
ary condition at the vortex core = 0, which physi- 



cally corresponds to a complete enclosure of the vortex 
flux within a finite-sized cylinder. Here we briefly dis- 
cuss other boundary conditions. For general value of 

the transition matrix elements (|45|) are modified at 

1 = (o-±l)/2: 



J 



U, 



(qi-l-rT,ki),(q 2 ,l,k 2 ) 



(qi,l-o-,ki),(q2,l,k 2 ) 



= cos(0) 
sin(0) 

= -sin(0) 
cos(0) 



2<rqiq 2 E x + E 2 Ca_fh_ 

7T Ei 2ir \k 2 



Tl-i 



Ei + E 2 , fh-k 2 
log UTfe 



al- 



4 V /.,M(.E_> - E t ) -C n ( ^ 



2aq 1 q 2 ( E x + E 2 

qiq 2 -p, + 2q\q 2 G 



\Jk\k 2 
* Ei + E 2 



\fk~xk 2 



e(a(k2-h)) 



Ei — Ei 



4 V '*iM(£j ■ E->) - C„ ( ^ 



9jl (hi 

2ir \k 2 
Ei + E 2 



al- 



l = 



(7+1 



Ei + E 2 , (ki-k 2 
log 



\/kik 2 



\Jk\k 2 



Q(a(ki-k 2 )) 



ki + k 2 



No new singularities have been introduced, and all re- 
sults from previous subsections remain qualitatively the 
same. However, there are quantitative changes to the 
vortex mass, finite temperature dissipation coefficient, 
etc. Therefore, the parameters that characterize vortex 
dynamics depend on the boundary conditions at the vor- 
tex core. We will not investigate further this dependence. 



APPENDIX F: REPRESENTATION OF THE 
FULL QUASIPARTICLE HAMILTONIAN 

Here we construct the representation of the Hamilto- 
nian (|33|1 in the basis of states JUJ, (|42|l . This is used to 
numerically diagonalize the Hamiltonian on a finite sam- 
ple with radius R. The basis wavefunctions (|4f|l and 
(I42|l are easily renormalized to unity on the finite sample 
if multiplied by yj n/R. With an appropriate hard-wall 
boundary condition momentum is quantized to zeros of 
the Bessel's functions: 









E\ = -( 


I- 1 - 




1 2R \ 


2 





fei?> 1 . (Fl) 



the section llVl we have (at the p = hkpx node): 

2ir R 

(Ei,h\H D \E 2 ,l 2 ) = [dtf, [dr (F2) 

sin< 







(v F cos 0)2 + (va sin0)2^>' 0)V ^ (r '' 
vfva {y/Ei)* ' \[E 2 



AR 



li'i'(ki,k 2 ) + 



</i^si»ii ( h - - j sign (h~^) -/;•./•/ 1' A- 1 . !»'•> I 



d<t> 



sin e 



i(l2—h)<t> 



(vf cos (f>) 2 + (va sin < 



^2 • 



where, like before, qi — ±f , Ei = qiki, and V and i" are 
indices of the Bessel's functions in the upper and lower 
spinor components respectively of the wavefunctions l|4f|) 
and l)42|) (which depend on !). Recall that v is gauge 
invariant, so that we cannot avoid its "elliptical" shape 
in the rescaled coordinate system. The integral over <j> 
can be calculated exactly; for I — l 2 — li > it is equal 
to: 



Note that the degeneracy with respect to angular momen- 
tum I is lifted in a finite system. In the chosen representa- 
tion, the full Hamiltonian l|33fl is a sum of a diagonal ma- 
trix that contains discretized energies E and the Dopplcr 
shift part whose matrix elements (Ei, li\H]j\E 2 , l 2 ) can 
be evaluated using the rescaled coordinate system intro- 
duced in section IV Al Assuming that the supercurrent v 
describes a 7r-flux localized at the origin, as discussed in 



f 



Vf + VA 



v F - va 



vf + va 



The new function 7 is: 



i l - {-if 



R 



vf > va 



l (f - (-1 )') , v F < V A 



7 h ,l 2 {ki,k 2 ) = / drJ h (kir)J h (k 2 r) . (F3) 
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Analytic expression for 7 is available only in the R — > 00 
limit (a special case of Weber-Schafheitlin formula) : 



r (jL±j2±l) 



k[ 2+1 r(j 2 +i)r(ii=£±l) 
fc H r(ii±^±l) 
4 1+1 r(ii+i)r(ia=^±i) 2 



2^1 



^1+^2 + 1 fa — + 1 . 



h + l; 



h+h + l h-h+l .J ri./fcj. 

9 ' 9 'I'M 



fel > fc 2 
fe 2 > fcl 



(F4) 



-^rsign (/1 + Z 2 + ±) sign (?i - Z 2 + ±) (-l)^ 211 
J_ cos (^^) log(fcR) 



h — h is an odd integer 
otherwise in the limit kR 3> 1 



It turns out that the formulas (|F1|) and (|F4|> pro- 
vide a very good approximation for our purposes. This 
approximation scheme uses the fact that for all states of 
interest kR 1 in the thermodynamic limit. The states 
with kR ~ 1 are not treated accurately when R is finite, 
but all such states collapse to an infinitely small range 
of energies when R — > 00. Then, the Hamiltonian repre- 
sentation can be constructed relatively fast. The results 
of diagonalization and vortex dynamics obtained within 
this approximation have been found virtually the same 
as the corresponding results with exact energy level dis- 
cretization and exact numerical calculation of 7 in finite 
systems. 



APPENDIX G: SEMI-ANALYTICAL 
CALCULATION OF TZ(Ae) 

Here we present a semi-analytical argument against 
dissipation and anomalous vortex dynamics caused by 
the Doppler shift. We make several assumptions about 
general properties of the exact wavefunctions ^ n (E, I) = 
(E,l\n) in the thermodynamic limit. First, we assume 
that the wavefunctions ^f n (E, I) are analytic functions of 
E everywhere except perhaps at E = g(e), and also at 
E = and/or e = 0. The function g(e) satisfies: 

5 (e) = -g(-e) (Gl) 
sign(c/(e)) = sign(e) , 

and gives the energy at which ty n (E,l) may be devel- 
oping a Dirac peak in the thermodynamic limit. These 
properties of g(e) are clearly visible from the numeri- 
cal spectra (the first one is protected by a symmetry). 
The numerics essentially suggests g(e) = e, and this is in 
agreement with scattering theory in the thermodynamic 
limit. No divergences occur at E — or e = 0, and these 
possible non-analyticities do not play an important role. 
The remaining assumptions follow from the requirement 
that the function IZ(Ae), which determines the vortex 
action, be finite for every finite Ae. All wavefunction 



singularities should be integrable; integrals of the wave- 
function (not its modulus squared) across the point of 
possible non-analyticity are convergent in the thermody- 
namic limit: 

g(e)+AE 

J dE \$ eti (E,l)\ < 00 . (G2) 

g(e)-AE 

This is certainly consistent with a possible Dirac peak. 
Also, certain sums over angular momentum channels 
should be either finite or lead to integrable singularities. 
All of these assumptions are based on the numerically 
obtained wavefunctions and spectra. 

We now analytically explore general properties of the 
function IZ(Ae) given by l|72() . and check whether any 
infra-red divergences might occur. A particular question 
is what the dependence on small Ae is. If TZ(Ae) ~ Ae^ 
with v < 0, then vortex dynamics at small frequencies 
may be anomalous and not simply governed by the vortex 
inertia. 

Let us introduce a convenient symbolic notation for the 
exact states \n) = |e, z), where e = e„ and i is a fictitious 
quantum number that replaces the angular momentum / 
of the states in absence of the Doppler shift. The purpose 
of this is only to explicitly reveal the role of energy e, 
which does not completely specify the exact states. Then 
we symbolically write (17211 as: 

A 

TZ(Ae) - E E / d£l / de 2 S(e 2 - £l - Ae)\U h2 \ 2 

nodes zi , t2_^ Q 
A A 

v x ,i = J2 J dE ^ J dE * ( G3 ) 

(Ei,li\ei,ii)(e2,i2\E2,h)r v (u>)U Elt i 1 -E 2 ,i 2 ■ 

Note that conventionally the summation over fictitious 
i\ and i 2 would be replaced by the appropriate density 
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of states factors. We integrate out ei in (|G3|I and write: 

Ac 



TZ(Ae) = E E / d£ l^ 



nodes ii ,*2 

A A 

u h2 -EE / dSi / dE2 (G4) 

Now we substitute individual terms from the expression 
(|45|l for transition matrix elements U e 1 ,i 1 ;E 2 .i 2 between 
the basis states, and check what happens as Ae — > 0. We 
start with potentially the most critical transitions, which 
involve the zero angular momentum channel. 



a. Transitions involving the zero angular momentum 

The transition matrix elements that we look first at 
are given in (|45|l for I = 



U Ei,l-<j;E 2 ,l <x - 



2a signal £2) Ex + E 2 
n Ex~E 2 



(G5) 



CV(<?i,g2) 

2vr 



E x 2 E x + E 2 , ( \E X 
:log 



E-, 



\/\E\E 2 



\Ei 



\E 2 



This expression, appropriate in the thermodynamic limit, 
appears formally infinite at E\ — E 2 , and we need to reg- 
ularize it on a finite sample in order to see the exact effect 
of this singularity. When the transition matrix elements 
U are calculated on a finite sample with radius R, i|G5|) 
is replaced by an oscillatory expression U^ R \E 2 — Ei), 
with period ~ R^ 1 and local average exactly given by 
(jG5fl . The discrete energies E (for \E\R^> 1) have spac- 
ing equal to the half-period of U^ R \ so that indeed in 
the R — > 00 limit the summation of over energies 
transforms into the integration over the "local average" 
(|H*5)l for every E 2 f E x . Exactly at E 2 = Ei = E it 
turns out that U^ R ' ~ \og{\E\R). However, this loga- 
rithmic divergence is innocuous, because its contribution 
to the sum over energies is _R _1 log(|-E|i?) — > in the 
R — ► 00 limit. In the following, we will symbolically ig- 
nore it in the thermodynamic limit by writing a principal 
part of (|G5|I . and focus on the divergent behavior close 
to Ei = E 2 . Note that there are no other divergences in 
this expression, in particular at E\ — > 0, E 2 — > or both. 

Consider U\. 2 in i|G4|) . We can change variables 
E\ = E — AE, E 2 = E. The transition matrix ele- 
ment <|G5|) diverges only at AE — > 0, so let us fix a 
small finite AE < Ae and imagine integrating out E 
first in (IG4II. For some value of E, one of the two wave- 



functions ^t-A£,i 1 (^^~ a ) an d ^ e ^ 2 {E 2 , 1) may become 
divergent, but not both - see Figure Q Based on numer- 
ics, we assume that a single wavefunction divergence is 
integrable, so that the integral over E is finite. As a func- 
tion of AE, this integral has a singularity inherited from 



the expression (|G5p at AE — 0, but its nature has not 
been changed and no other singularity at any AE < Ae 
has been introduced. This is important, because if we 
now proceed with integrating out AE, the singularity at 
AE = behaves exactly the same way as if we were sim- 
ply integrating l|G5() . However, it is not difficult to check 
that such an integral: 



A-Eq 



dA£?( const. 



AE 



const. log(A£ 2 ) 



(G6) 



-AE 



evaluates to a finite value (P denotes the principal part 
due to finite sample regularization) . 

What remains to be examined is how sending Ae to 
zero may affect this picture. Generally, the integra- 
tion over E discussed above picks divergences of both 
wavefunctions in l|G4l) . so that the resulting function 
of AE may be non-analytic or even divergent at some 
AE ~ Ae. In particular, if the only divergences of the 
wavefunctions are Dirac peaks (as it seems), then the 
resulting function of AE is divergent as a Dirac peak at 
AE = g(e)—g(e—Ae) ~ Ae; other finite non-analyticities 
in the wavefunctions do not matter. Now, integrating 
out AE is sensitive to such a divergence and the final re- 
sult (which should be finite for Ae ^ 0) may end up 
scaling as U(AE ~ Ae) ~ ^ for small Ae. How- 
ever, this can actually happen only if we introduce some 
non-zero chemical potential (i. Otherwise, we always ap- 
proach the Ae ~ E 2 — E\ — * limit only together with 
e ~ E\ + E 2 — > 2/i = 0, and no divergence emerges from 
(|G4p (for details, see analytical calculations in absence 
of the Doppler shift - this is equivalent to reducing the 
present exact wavefunctions to pure Dirac peaks, and 
such divergences do not give rise to dissipation without 
a non-zero chemical potential). 

In conclusion, the transitions that involve the zero an- 
gular momentum, which are the most critical for appear- 
ance of anomalous vortex dynamics, seem to qualitatively 
follow the same behavior as if there were no Doppler shift. 



E, =E-AE 



E 2 =E 



AE 



~Ae 



g(e-Ae) 



g(e) 



FIG. 7: Integration of (7i,2 in l|Q4(l over E. Ticks mark 
the points at which the wavefunctions are possibly divergent 
(K-Ae, H {Ei,l - cr) = *:_ Ae>il (B - AE,l - a) at the left 
tick, and * £ ,i 2 (_E2,0 = >I y e ,i 2 (i?, I) at the right tick). As we 
integrate over E, for fixed A_E, the dashed lines move from 
left to right in this diagram, and cross over the singularities 
in the wavefunctions. Only for one value of AE ~ Ae both 
singularities can be crossed at the same time. 
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Transitions that do not involve the zero angular 
momentum 

Consider first the "most divergent" terms in (|45|l for 

+J 

2 - 

<j + r 



U El ,i-a;E 2 ,l oc 4^1^1^21-5(^2 - £i)sign 

' (G7) 

and substitute them into I|(t4(1 . The sign function takes 
care of the signs for positive or negative I in Ij45(l . Ig- 
noring r v (a;) and the vector structure of U, we find: 

ft,2« E E / d ^-A Ml (^-^R >i2 (£,o 



|£|sign i 



<r + l 



(G8) 



According to our assumptions from the beginning of this 
appendix, the energies E at which the wavefunctions 
\J/ e i(i?, I) are non-analytic do not depend on angular mo- 
mentum I. Therefore, we can first sum up I - the result 
should be finite, and have the same analytical properties 
as the original expression. Then, we integrate out E. 
As long as Ae ^ 0, the two wavefunctions in the integral 
cannot diverge at the same value of E, so that, according 
to our assumptions, the integral over E is finite. Exactly 
for Ae = 0, the result may be infinite, and it is S(Ae) if 
the wavefunction divergence at E ~ e is a Dirac peak. 



This is essentially the same situation as in the case with- 
out the Doppler shift: <5(Ae) is killed by the factors of 
Ae in the vortex action. 

Now consider the remaining portion of l|45l) for I ^ 



^-ji and substitute: 



U E u l-a;E 2 ,l OC | — r 



Ei+E 2 



( sign ( I — 



^/\E 1 E 2 
a + l 



a(|£ 2 | - \E X \) 



into Ui.2 from (|G4J| . The sign function selects whether 
the argument of the step O-function should be a(k 2 — ki) 
for positive I or a (hi — k 2 ) for negative I (see l|45[) . and 
recall k = \E\,q = sign(£')). We could roughly follow 
the previous argument, although the summation over I 
is not entirely trivial - it may introduce a singularity 
at E\ — E 2 . Since the wavefunction amplitudes cannot 
grow indefinitely with I, this new singular behavior is in 
the worst case: 



Ei + E 2 
Ei — E 2 



(G9) 



(see for example l|49[l). After summation over I we are 
left with an integral over Ei and E 2 , which has a very 
similar structure to the one that we considered in the pre- 
vious subsection. By repeating the same argument, we 
arrive to the conclusion that we could obtain anomalous 
vortex dynamics only with a non-zero chemical potential. 
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